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ABSTRACT 


Tilmvise  condensation  measurements  of  steam  were  made  on 
horizontal  finned  tubes  under  vacuum  and  near-atmospheric 
conditions.  Data  were  obtained  for  copper  tubes  with  fins 
of  rectangular,  triangular,  trapezoidal,  and  parabolic  cross 
sections,  and  for  a  commercially-available  finned  tube.  A 
stainless  steel  finned  tube  was  also  tested  to  investigate 
the  effect  of  thermal  conductivity. 

jlaximum  enhancements  of  about  4.8  were  obtained  under 
vacuum  conditions,  and  about  6.9  at  atmospheric  pressure, 
compared  to  a  smooth  tube  having  an  outside  diameter  equal 
to  the  root  diameter  of  the  finned  tubes.  The  optimum  fin 
spacing  was  found  to  be  about  2.0  mm  for  rectangularly 
shaped- fins  with  a  fin  thickness  of  1.0  mm,  and  fin  height 
of  0.5  and  1.5  mm.  Fins  with  a  parabolic  shape  were  shown 
to  perform  better  than  fins  of  rectangular  shape,  and  fins 
were  shown  to  degrade  the  performance  of  stainless  steel 
tubes.  The  effects  of  vapor  shear  were  shown  to  have  only  a 
small  influence  on  the  steam-side  heat-transfer  coefficient. 
A  theoretical  model  proposed  by  Rebb  et  al.  [25]  was  found 
to  underpredict  the  experimental  data.  Several  suggestions 
to  modify  this  model  are  described.  r  '  -•  '• 
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I.  IHTRODOCTIOM 


A.  BACKGBOUBD 

As  costs  for  material  and  energy  continue  to  rise,  the 
need  for  smaller,  more  efficient  heat  exchangers  continues 
to  grow.  For  those  applications  involving  marine  vehicles, 
size  and  weight  limitations  dictate  the  use  of  small, 
highly-efficient  heat  exchangers.  Reducing  the  size  and 
weight  of  steam  condensers  used  aboard  a.  S.  Navy  ships  for 
propulsion  and  electrical  power  generation  would  result  in 
lower  material  costs,  and  help  to  aleviate  the  cramped 
conditions  so  typical  of  machinery  spaces.  Thus,  there  is  a 
strong  motivation  for  continued  study  in  this  area. 

The  effectiveness  of  condensers  is  limited  by  the  water¬ 
side,  vapor-side,  and  wall  thermal  resistances  of  the 
condenser  tubes.  Reducing  any  one  of  these  will  contribute 
to  increased  heat- transfer  performance,  and  smaller  physical 
size  of  condensers.  Methods  of  enhancing  the  vapor-side 
coefficient  include  the  use  of  "roped"  tubes,  fluted  tubes, 
drainage  strips  attached  on  the  tabes,  finned  tubes,  and 
coatings  applied  to  enhance  dropwise  condensation.  This 
thesis  concentrates  soley  on  finned  tubes. 

Since  the  late  1940s,  externally-finned  tubes  have  been 
used  to  increase  the  vapor-side  heat-transfer  coefficient  of 
tubes  used  in  refrigeration  systems;  but  condensers  used  in 
steam  systems,  such  as  shipboard  propulsion  plants,  continue 
to  use  smooth  tubes.  The  high  surface  tension  of  water, 
which  leads  to  its  tendency  to  flood  the  area  between  fins, 
has  resulted  in  a  widely-held  belief  that  such  tubes  are 
inappropriate  for  use  in  steam  systems.  Recent  studies 
[1,2],  however,  have  shown  that  finned  tubes  can  signifi¬ 
cantly  enhance  the  heat-transfer  rates  in  such  systems. 
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The  theoretical  treatment  of  condensation  on  finned 
tubes  is  extremely  complex  due  to  the  large  number  of  vari¬ 
ables  and  physical  mechanisms  involved.  The  interaction  of 
gravitational  and  surface-tension  forces  lead  to  complex 
three-dimensional  flow  patterns,  which  are  further  dependent 
on  fin  spacing,  height,  and  thickness.  Other  variables 
include  heat  flux,  vapor  shear,  tube  diameter,  fin  shape  and 
fluid  properties  just  to  name  a  few.  In  view  of  the  above, 
any  theoretical  models  will  require  numerous  simplifying 
assumptions,  and  require  complex  computer  solutions 
involving  implicit  numerical  techniques  [3].  To  confirm  the 
validity  of  theoretical  models,  reliable  experimental  data 
which  cover  a  wide  range  of  relevant  parameters  must  be 
obtained.  The  availability  of  such  data  may  then  lead  to 
the  development  of  fairly  simple  experimental  correlations, 
which  could  be  used  in  the  design  stage  to  predict  the  heat- 
transfer  performance  of  finned  tubes. 

This  thesis  effort  is  a  continuation  of  research  being 
conducted  at  the  Naval  Postgraduate  School  (NPS)  under  a 
grant  from  the  Rational  Science  Foundation.  The  basic  test 
apparatus  used  to  ccUect  experimental  data  was  built  by 
Krohn  [4].  Graber  [5]  provided  the  majority  of  instrumenta¬ 
tion,  and  took  preliminary  data  as  the  system  experienced 
problems  with  non-condensing  gases  and  partial  dropwise 
condensation  on  copper  tubes.  Poole  [6]  made  further 
improvements  on  the  apparatus  as  well  as  on  the  instrumenta¬ 
tion  and,  most  importantly,  assured  a  leak-free  apparatus. 
Unfortunately,  he  did  not  have  sufficient  time  to  produce 
useful  data,  mostly  due  to  the  considerable  time  spent  in 
systematically  locating  and  fixing  leaks,  and  due  to  the 
partial  dropwise  condensation  problem  that  had  not  been 
solved.  Using  this  system,  Georgiadis  [7]  was  finally  able 
to  obtain  complete  filmwise  condensation  on  copper  tubes. 
He  took  data  on  a  number  of  finned  tubes  with  fins  of 
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rectangular  shape,  as  well  as  on  saiooth  tubes.  The  repeat¬ 
ability  of  data  obtained  by  Georgiadis  demonstrated  the 
accuracy  of  the  test  apparatus  and  associated  instrumenta¬ 
tion  which  was  used  essentially  without  modification  for 
this  investigation. 

The  overall  objectives  of  the  present  program  at  NPS 
includes  the  testing  of:  (a)  tubes  with  rectangularly- 

shaped  fins  to  find  the  best  fin  spacing,  thickness,  and 
height,  (h)  tubes  with  various  fin  shapes  to  find  the  best 
geometry  that  will  maximize  heat  transfer,  (c)  tubes  with 
different  tube-metal  thermal  conductivity,  and  (d)  the 
effect  of  vapor  shear  on  finned  tubes.  Georgiadis  [7] 
tested  a  total  of  25  tubes,  which  were  divided  into  five  fin 
spacings  (0.5,  1.0,  1.5,  2.0  and  4.0  mm),  five  fin  thick¬ 
nesses  (0.5,  0.75,  1.0,  1.25,  and  1.5  mm),  and  two  fin 
heights  (1.0  and  2.0  mm).  At  the  conclusion  of  Georgiadis* 
test  program,  eight  new  tubes  with  fin  heights  of  0.5  and 
1.5  mm,  fin  spacings  of  1.0,  1.5,  2.0,  and  4.0  mm,  and  a 
fin  thickness  of  1.0  mm,  remained  to  be  tested  to  complete 
the  sequence  of  tubes  with  rectangularly- shaped  fins.  The 
testing  of  these  tubes  and  additional  tubes  toward  objec¬ 
tives  (b)  ,  (c)  ,  and  (d)  were  the  primary  goals  of  this 

thesis  effort  as  listed  in  the  next  section. 

B.  OBJECTIVES 

The  main  objectives  of  this  thesis  are  as  follow: 

1.  Take  data  on  tubes  with  rectangularly-shaped  fins  of 
various  fin  spacings  and  fin  heights  to  augment 
previous  data  [7  ], 

2.  Take  data  on  tubes  with  fins  of  different  shapes 
(triangular,  trapezoidal,  parabolic,  etc.). 
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ake  data  on  commercially-available  finned  tubes. 


Take  data  on  the  "optiBum rectangularly-shaped 
finned  tube  at  different  vapor  velocities. 

Take  data  on  a  stainless  steel  tube  with  "optimum," 
rectangularly -shaped  fins,  and 

Develop  a  preliminary  correlation  based  on  data  for 
tubes  with  rectangularly-shaped  fins. 


II.  THEORETICAL  TREATHEHT  OF  COMDENSATIOM  ON  BORIZONTAL 

115 NED  TUBES 


A.  FILfl  CONDENSATION 

Rhen  filmwise  condensation  of  a  vapor  takes  place  on 
smooth  tubes,  a  thin  layer  of  condensate  is  formed  vhich 
thickens  with  increasing  distance  around  the  perimeter  of 
the  tube.  This  condensate  layer  creates  a  thermal  resis¬ 
tance,  which  can  limit  the  heat- transf er  performance  of  the 
tube.  This  film  thickness  and  its  thermal  resistance  can  be 
reduced  by  the  use  of  external,  radial  fins  which,  in  addi¬ 
tion  to  a  surface  area  increase,  promote  surface-tension 
effects. 

In  1984,  Yau  et  al.  [1]  measured  the  enhancement 
provided  by  copper  finned  tubes  over  smooth  tubes  for  film- 
wise  condensation  of  steam.  Similar  experiments  by 
Wanniarachchi  et  al.  [2]  also  in  1984  confirmed  that  the 
observed  enhancements  were  greater  than  could  be  explained 
by  the  increased  surface  area  alone.  This  additional 
enhancement  may  be  a  result  of  the  surface-tension  forces 
which  act  to  thin  the  condensate  film.  Figure  2.1  schemati¬ 
cally  depicts  this  phenomenon. 

The  effect  of  surface  tension  on  pressure  at  the  inter¬ 
face  between  a  liguid  and  vapor  is  inversely  proportional  to 
the  radius  of  curvature  of  the  interface.  If  a  surface  has 
two  radii  of  curvature  at  right  angles  (i.e. ,  r^  and  r  2)  r  it 
can  be  shown  that: 


AT 


(2.1) 
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o  =  surface  tension  of  condensate,  and 
A  F  =  the  pressure  difference. 


Tin  Condensate 


Figure  2.1  Schenatic  of  Condensate  Profile  at  Upper 

Tube  Surface. 

For  the  case  of  a  finned  tube,  the  radius  of  curvature 
around  the  outer  fin  perimeter  is  very  large  compared  to  the 
radius  of  curvature  around  the  fin  cross-section  profile 
(Figure  2.1).  Peferring  back  to  equation  (2.1)  ,  it  can  be 
seen  that  the  smaller  radius  of  curvature  term  will  domi¬ 
nate,  so  surface-tension  effects  around  the  fin  perimeter 
may  be  neglected. 

Because  of  the  convex  shape  of  the  condensate  film  at 
point  A,  the  pressure  within  the  film  at  this  point  is 
greater  than  the  surrounding  vapor  pressure.  In  a  similar 
manner,  the  pressure  at  point  B  is  less  than  the  surrounding 
vapor  pressure  owing  to  the  concave  shape  of  the  condensate 
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film.  The  relatively  flat  shape  of  the  condensate  film  at 
point  C  leads  to  a  pressure  essentially  euual  to  the 
surrounding  vapor  pressure.  These  pressures  are  given  by: 


P.  =  P  +  £_ 
A  V 


(2.2) 


P 


p. 


(2.3) 


r 


r. 


r 


V 


(2.4) 


where 


surrounding  vapor  pressure, 

liquid  pressure  at  points  A,  E,  and  C,  and 


Cbt  Tc 


radius  of  curvature  of  the  condensate 
film  at  points  A,  B,  and  C. 


At  point  R,  the  radius  of  curvature  is  small,  so  the  pres¬ 
sure  at  point  D  is  less  than  the  pressure  at  point  C. 
Further,  from  equations  (2.2)  and  (2.3),  the  pressure  at 
point  A  is  greater  than  at  point  B.  In  reality,  the  pres¬ 
sure  gradient  within  the  condensate  film  varies  along  the 
hnight  of  the  fin  due  to  the  continuously  varying  radius  of 
curvature  from  the  fin  tip  down  to  the  base  [3];  however,  to 
simplify  tl.e  treatment  of  condensate  flow  toward  the  fin 
base,  tlie  pressure  differences  between  points  A  and  B  and 
points  C  and  B  may  be  simply  written  as: 
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vhere  j 

ap^B'  ~  pressure  difference  between  points  A  and  B, 

and  points  C  and  B. 


As  can  be  seen,  these  pressure  differences  are  positive, 
rts’jitin<j  in  condensate  flow  toward  point  B.  Due  to  the 
relatively  larye  mass  of  condensate  at  point  B,  gravita¬ 
tional  forces  dominate  over  surface-tension  forces  causing  a 
flow  of  condensate  around  the  tube  perimeter  at  the  fin 
base.  In  this  manner,  a  condensate  run-off  channel  is 
formed  at  point  B,  and  the  improved  drainage  of  condensate 
thins  the  film  between  fins  and  on  the  fin  surface.  This 
thinning,  in  turn,  reduces  .  the  thermal  resistance  through 
the  condensate  film,  thus  producing  an  enhancsaent  in  addi- 
tioii  to  the  gain  in  surface  area. 

It  sl.oull  te  stressed  at  this  point  that  the  shape  of 
the  condensate  film  is  highly  dependent  on  the  fin  geometry. 
Clearly,  tJjen,  the  fin  shape  must  be  taken  into  considera¬ 
tion  in  order  to  maximize  the  beneficial  effects  of  surface 
tension. 

Tlie  above  gains  may  be  partially  or  totally  offset, 
however,  Ly  the  tendency  of  condensate  to  flood  the  area 
between  fins,  especiallv  on  the  xower  part  of  the  tubes. 


The  extent  to  which  the  film  floods  the  tube  is  defined  by 
the  condensate  retention  angle  (0),  as  shown  in  Figure  2.2. 
In  the  flooded  portion  of  the  tube,  the  relatively  thick 
condensate  film  increases  the  thermal  resistance,  thereby 
leading  to  a  degradation  in  the  heat-transf er  performance. 
The  condensate  retention  angle  is  highly  dependent  on  fin 
spacing,  so  very  closely-spaced  fins  may  lead  to  completely 
flooded  tubes. 

A  number  of  studies  have  been  conducted  to  investigate 
the  performance  of  finned  tubes,  and  these  are  presented  in 
the  next  section. 

B.  COHDElfSATE  RETEHTION 

In  1946,  the  first  measurements  of  condensate  retention 
were  made  by  Katz  et  al.  [8].  These  measurements  were  made 
under  static  conditions  (i.e.,  no  condensation  taking 
place),  using  water,  aniline,  acetone,  and  carbon  tetrachlo¬ 
ride  on  a  number  of  tubes  with  different  fin  densities  (276 
to  984  fins/m),  and  fin  heights  (1,2  to  5.7  mm).  It  was 
shown  that  as  much  as  1005?  of  the  tube  surface  could  be 
flooded  by  retained  condensate,  depending  mainly  on  the 
ratio  of  surface  tension  to  liquid  density  and  on  the  fin 
spacing. 

Eight  years  later,  Gregorig  [9]  recognized  that  surface- 
tension  effects  could  play  an  important  role  on  a  vertical, 
fluted  surface  which  uses  a  minimum  radius  of  curvature  at 
the  flute  tip,  which  gradually  increases  toward  the  trough. 
As  discussed  in  section  A  above,  this  variable  radius  of 
curvature  results  in  a  condensate  film  pressure  that 
decreases  toward  the  trough.  In  this  manner,  the  condensate 
layer  on  the  convex  region  of  the  flute  is  thinned  signifi¬ 
cantly,  thus  pushing  the  condensate  into  the  troughs  (Figure 
2.3)  where  gravitational  forces  result  in  enhanced  drainage. 


External  Diameter  of  fins 


Figure  2.2  Scheaatic  of  Condensate  Retention 
on  Finned  Tubes. 
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Since  the  heat- tL'ansfer  coefficient  is  inversely  propor¬ 
tional  to  the  thickness  of  the  condensate  layer,  heat- 
tiausfer  performance  at  the  trough  is  degraded,  but  to  a 
lesser  degree  than  the  enhanceoent  experienced  by  the  crest. 
The  resulting  average  heat-transfer  coefficient,  therefore, 
shows  a  significant  enhancaoent. 


THIN  riLM  MEGION 


Figure  2.3  Cross  Section  of  the  Gregorig  Surface 


In  1981,  Rudy  and  Webb  [10]  made  measurements  of  conden¬ 
sate  retention  angles  on  finned  tubes  with  three  different 
fin  densities  (748,  1024,  1378  fins/neter)  using  water, 
P.-11,  and  n-pentane  as  the  working  fluid.  Later,  in  1983, 
Pudy  and  Webb  [11]  developed  an  analytical  model  to  predict 
tha  fraction  of  tube  surface  that  is  flooded  during  conden¬ 
sation  on  a  horizontal,  integral-fin  tube.  They  found  that 
the  vertical-rise  height  of  condensate  on  a  finned  tube  was 
the  same  as  that  obtained  on  a  vertically  oriented,  flat 
finned  plate  that  was  obtained  by  splitting  and  unrolling  an 
identical  finned  tube.  Based  on  this  observation,  they  used 
capillary  eguations  that  predict  the  liquid  rise  on  a 
vertical  U-shaped  channel,  and  assuming  negligible  vapor 
shear,  they  predicted  the  condensate  retention  angle  to  be: 
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e  =  fin  height, 
t  =  fin  thickness, 

Pj  =  density  of  condensate, 
g  =  acceleration  of  gravity, 
s  =  fin  spacing,  and 
Dq  =  root  diameter  of  tube. 

Tills  result  shows  that  the  condensate  retention  angle 
increases  both  with  increasing  fin  density,  and  with 
increasing  surface  tension-to-density  ratio.  It  should  be 
kept  in  mind,  however,  that  the  model  was  based  on  a 
vertical  surface,  so  it  should  only  be  used  for  angles  less 
than  about  30  degrees.  For  angles  below  this,  experimental 
results  involving  the  use  of  water,  E-ll,  F-12,  ammonia,  and 
n-pentane  were  predicted  to  within  10  percent.  It  is  worth 
mentioning  that  equation  (2.7)  had  been  derived  in  1982, 
according  to  Russian  literature,  by  Rifert  [  12]. 

Like  Rudy  and  Hebb  [10],  Owen  et  al.  [13]  also  recog¬ 
nized  the  need  to  consider  condensate  retention  while 
analyzing  finned  tubes.  In  order  to  correct  for  this,  an 
assumption  was  made  that  the  condensate  retention  angle  was 
independent  of  condensation  rate,  so  a  static  analysis  was 
performed.  k  simple  force  balance  between  surface  tension 


ar. J  gravitational  forces  resulted  in  an  equation  for  the 
condensate  retention  angle  as  shown  below: 
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=  overall  diameter  of  fins. 


This  equation  is  the  sane  as  equation  (2.7),  except  that 
equation  (2.8)  is  independent  of  fin  thickness  (t)  . 

In  1983,  Honda  et  al.  [3]  performed  a  theoretical  anal¬ 
ysis  to  determine  the  condensate  retention  angle.  Using  an 
iterative  numerical  scheme,  they  found  the  solution  to  agree 
witli  equation  (2.8)  obtained  earlier  by  Owen  et  al. .  Using 
data  of  their  own,  Honda  et  al.  verified  a  close  agreement 
between  the  predicted  and  experimental  values  of  the  conden¬ 
sate  retention  angle. 

Yau  et  al.  [1]  conducted  experiments  to  determine  the 
effects  of  removing  retained  condensate  from  finned  tubes  by 
installing  thin  metal  drainage  strips  attached  edgewise  to 
the  bottom  tube  surface.  Experiments  were  conducted  where 
the  vapor  velocity  and  fin  pitch  were  varied  on  finned  tubes 
with  and  without  drainage  strips.  They  showed  a  significant 
reduction  in  the  in  the  condensate  retention  angle  when  the 
finned  tubes  were  fitted  with  drainage  strips.  Condensate 
retention  angles  for  condensation  of  steam,  ethylene  glycol, 
and  1.-113  on  finned  tubes  with  drainage  strips  were  found  to 
fit  the  empirical  relation  listed  below  to  within  ±  10 
percent: 


>).  =  cos 
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In  1985,  Rudy  and  Webb  [14]  expanded  their  1983  nodel  in 
order  to  predict  the  condensate  retention  angle  for  hori¬ 
zontal  finned  tubes  with  fins  of  arbitrary  shape.  As 
before,  their  model  was  based  on  capillary  equations  that 
predict  the  amount  of  liquid  rise.,  and  negligible  vapor 
shear.  The  resulting  equation  for  the  condensate  retention 
angle  is  given  by; 
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vhere 


Pj^  =  wetted  perimeter  of  fin  cross  section, 
t^  =  fin  base  thickness,  and 

Ap  =  profile  area  of  fin  over  fin  cross  section. 


In  order  to  test  the  validity  of  equation  (2.10),  it  was 
compared  to  experimental  data  for  horizontal  finned  tubes 
with  fins  of  primarily  trapezoidal  cross  section,  and  fin 
densities  ranging  from  630  to  1614  fins/m.  Steam,  R-11,  and 
n-pentane  were  used  as  the  condensing  fluids.  Equation 
(2.10)  was  shown  to  predict  the  experimental  data  to  within 
±  ID  percent.  Note  that  equation  (2.10)  reduces  to  equation 
(2.8)  for  fins  of  rectangular  cross  section. 
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C.  THEORETICAL  MODELS 


In  1947,  Beatty  and  Katz  [15]  performed  a  number  of 
experiments  where  various  refrigerants  were  condensed  on 
single  finned  tubes.  To  predict  the  film  coefficients,  they 
started  with  the  Nusselt  eguations  for  condensation  on  a 
horizontal  tube  and  on  a  vertical  plate,  and  with  the 
respective  proportions  of  horizontal  tube  area  (between 
fins)  and  vertical  fin  area,  they  calculated  the  total  heat 
transfer  rate.  From  this,  an  average  heat-transfer  coeffi¬ 
cient  was  obtained  based  on  an  eguivalent  tube  diameter. 
Their  final  expression  is  simply  the  Nusselt  equation  for  a 
smooth  horizontal  tube,  but  with  the  tube  diameter  replaced 
with  the  eguivalent  tube  diameter.  The  leading  coefficient 
was  modified  to  fit  their  experimental  data  as  shown  below: 
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where 


Aq  =  surface  area  of  smooth  tube, 
h  =  average  vapor-side  heat-transfer  coefficient, 

30 


O  V  •  * 


Kf  =  thermal  conductivity  of  condensate. 


=  density  of  vapor, 

hfg  =  specific  enthalpy  of  vaporization, 

=  viscosity  of  condensate, 

AT  =  vapor- side  temperature  drop, 

Dg  =  equivalent  tube  diameter, 
n  =  fin  efficiency, 

A„  =  surface  area  of  smooth  tube, 

s 

Aj  =  total  surface  area  of  finned  tube,  and 
A  =  effective  area  of  finned  tube. 

Accuracy  of  this  equation  was  claimed  to  be  better  than  11 
percent  for  a  wide  variety  of  nonaqueous  fluids.  It  should 
be  noted,  however,  that  this  model  ignores  the  effect  of 
surface-tension  forces;  so,  it  is  valid  only  for  low- 
surface-tension  fluids  such  as  refrigerants. 

Analytical  and  experimental  studies  of  condensation  on 
horizontal  tubes  with  tra pezoidally-shapei  fins  were 
performed  by  Zozulya,  Karkhu,  and  Borovkov  [16,17]  some 
years  later  in  the  1970s.  Their  experiments  confirmed  the 
need  to  consider  surface-tension  forces  in  addition  to 
gravitational  forces  when  developing  mathematical  models. 
The  analytical  solutions  were  based  on  the  following  assump¬ 
tions:  1)  the  thin  condensate  film  on  the  fins  was  treated 

as  a  laminar  boundary  layer  with  a  pressure  gradient  along 
the  fin  profile  caused  by  surface-tension  forces;  2)  the 
effect  of  gravitational  and  inertial  forces  on  the  motion  of 
the  film  along  the  side  surfaces  of  the  fins  into  the 


condensate-fiLled  trough  was  neglected;  3)  the  motion  of 
condensate  in  the  trough  area  is  laminar  and  produced  by 
gravity;  and  4)  no  condensation  takes  place  on  the  flooded 
portion  of  the  tube.  Differential  equations  were  obtciined 
for  the  height  of  condensate  between  fin  bases,  and  tempera¬ 
ture  distribution  along  the  fin  height.  With  the  aid  of  a 
computer,  solutions  for  these  equations  were  obtained  using 
numerical  methods,  leading  to  the  expressions  below: 
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where 


=  condensate  flow  rate, 

=  fin  tip  width, 

=  fin  semivertex  angle,  and 

=  dimensionless  depth  of  condensate 
between  fins. 


Although  the  above  equations  appear  to  be  fairly  simple,  it 
must  be  kept  in  mind  that  rather  complex  computer  solutions 


are  required  to  obtain  values  for  Z.  To  check  the  validity 
of  equation  (2.15),  (2.16)  and  (2.17),  experiments  were 
performed  with  four  different  finned  tubes  to  condense  both 
steam  and  R-113  vapor.  Agreement  to  within  5  percent  was 
reported. 

In  1977,  Nader  [18]  presented  an  analytical  model  to 
predict  condensation  on  a  vertical,  single  fin  attached  to  a 
horizontal  tube.  This  model  differed  from  those  discussed 
above  in  that  the  fin  temperature  was  allowed  to  vary  along 
the  height  of  the  fin.  For  the  special  case  where  the  fin 
temperature  is  constant,  the  Nusselt  equation  was  obtained. 
A  "condensation  efficiency"  was  then  proposed  to  account  for 
a  variable  fin  temperature,  and  was  defined  as: 
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where 


=  constant  of  proportionality,  and 

k  =  thermal  conductivity  of  fin. 
m  •* 


TTith  this  definition,  the  rate  of  condensate  formation  and 
the  rate  of  heat  transfer  (with  variable  fin  temperature) 
were  found  to  be: 
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Gj  =  rate  of  condensate  formation,  and 
Q  =  total  heat-transfer  rate. 


Thus,  with  a  table  of  values  for  and  obtained  numeri¬ 
cally,  the  rate  of  heat  transfer  and  condensate  formation 
could  be  calculated  by  a  fairly  simple  expression  provided 
the  condensation  efficiency  is  greater  than  0.8.  However, 
for  values  below  0.8,  the  accuracy  of  equations  (2.21)  and 
(2.22)  diminishes  quite  rapidly. 

In  1979,  Patankar  and  Sparrow  [19]  treated  the  case  of  a 
fin  attached  to  a  cooled  vertical  plate  or  a  cylinder  as  a 
three-dimensional  problem.  A  thin  fin  was  assumed,  and 
temperature  variations  acrass  the  thickness  were  ignored  as 
was  lengthwise  conduction.  Additionally,  all  heat  transfer 
was  assumed  due  to  condensation  only.  Compared  to  the  ideal 
fin  solution,  results  obtained  by  Patankar  and  Sparrow  indi¬ 
cated  that  a  significant  error  existed  when  using  the  ideal 
fin  solution  on  all  but  very  long  fins.  Based  on  their 
results,  the  following  relations  were  proposed: 


q  =  I.Oft  k  ^  '-f  '‘f 

™  e  bfj,  <>  Pf 


n  =  n.fti7i  ^’fr  [*'  ^  ^  ^ 


•'m  *■ 

7  kf  e2 


(2.  23) 


2.3-  -.^''Sr  n^''2 

i-1.  ''  "f  '=»  M2.  24) 

r  bfj,  p  pf^  2  kf 


where 


g  =  heat  flux,  and 


2  =  axial  coordinate. 


According  to  Patankar  and  Sparrow,  these  equations  should  be 
valid  for  most  practical  applications;  however,  no  experi¬ 
mental  data  were  available  at  the  time  to  confirm  this. 

In  1980,  Rifert  [20]  analyzed  condensation  of  vapor  on 
horizontal  finned  tubes  enhanced  by  the  effect  of  surface- 
tension  forces.  In  his  analysis,  he  divided  the  tube  into 
flooded  and  unflooded  zones,  and  solved  a  two-dimensional 
form  of  the  energy  equation  for  each  zone.  The  mean  heat 
flux  was  then  determined  by  integrating  over  each  zone  and 
the  tube  length.  In  cases  where  condensate  is  retained  in 
more  than  half  of  the  tube  perimeter,  Rifert  points  out  that 
a  three-dimensional  form  of  the  energy  equation  must  be 
used.  Solutions  to  these  equations  revealed  that,  in  most 
cases,  the  fin  temperature  is  very  nonuniform,  so  the  mean 
integral  heat-transfer  coefficient  and  its  correlations 
should  not  be  used.  The  above-mentioned  zone-by-zone  anal¬ 
ysis  was,  therefore,  recommended. 

Based  on  their  study  of  condensate  retention  mentioned 
earlier  in  section  B,  Rudy  and  Webb  [10]  proposed  that  the 


Beatty  and  Kate  [15]  model  be  modified  to  account  for 
condensate  retention,  as  shown  below: 


>1 


(2.  25) 


where  h_-,  is  computed  using  eguation  (2.11).  This  equation 
neglects  any  heat  transfer  through  the  flooded  portion  of 
the  tube.  As  a  result,  it  was  shown  to  underpredict  the 
average  heat-transfer  coefficient  of  condensing  R-11  by  as 
much  as  30  percent  when  a  significant  amount  of  retained 
condensate  was  present. 

Using  the  condensate  retention  angle  they  had  developed 
earlier  (see  section  B)  ,  Owen  et  al.  [13]  divided  a  hori¬ 
zontal  finned  tube  into  an  upper  unflooded  section,  cind  a 
lower  flooded  section.  Unlike  Rudy  and  Pebb  [10],  however, 
they  accounted  for  heat-transfer  through  the  flooded  portion 
of  the  tube  by  noting  that  condensation  occurs  on  both  the 
surface  of  the  retained  condensate  and  the  fin  tips.  An 
effective  thermal  resistance  was  obtained  over  this  portion 
of  the  tube  by  assuming  parallel  paths  for  heat- transfer 
through  the  fins  and  retained  condensate.  The  respective 
thermal  resistances,  which  were  added  in  parallel,  were  then 
added  in  series  with  the  tube  wall  thermal  resistance.  The 
resulting  heat-transfer  coefficient  was  then  found  to  be: 


h 


(2.  26) 
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vhere  is  computed  using  equation  (2.11)  ,  and 


k 


+  1 


Wf  _ 

kf^  Pf  rpf  -  pv>  n  >ifp. 


1/4-1  “1 


)  ] 


(2. 27) 


k  =  effective  thermal  conductivity 
erf 


This  equation  was  found  to  precict  all  the  available  data 
with  an  accuracy  of  better  than  30  percent. 

Pudy  and  Webb  [21]  used  equation  (2.7)  a  short  time 
later  while  developing  an  equation  for  the  heat-transfer 
coefficient.  The  model  for  this  undertaking  was  based  on 
surface-tension- driven  radial  flow  of  condensate  on  the  fin 
surface,  with  a  linear  pressure  gradient  along  the  fin 
height,  and  gravity-drained  flow  of  condensate  in  the 
channel  between  fins.  The  Nusselt  equation  for  horizontal 
tubes  was  used  for  the  tube  area  between  fins,  while  the  fin 
surface  was  treated  by  an  equivalent  gravity  model  developed 
by  Rebb  et  al.  [22]  and  Rudy  [23]  earlier.  As  before,  vapor 
shear  and  heat  transfer  through  the  flooded  portion  were 
neglected,  yielding  a  heat-transfer  coefficient  of: 


n.7'’5 


0.043  n 


3  2 

kf  Pf  iifp  n 
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pf  rj,  AT  “ 
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(2.  28) 
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1  =  length  of  tube, 

=  surface  area  of  tube  between  fins,  and 
Ajt  =  fin  surface  area. 


This  expression  provided  an  accuracy  of  better  than  1055  for 
condensation  of  E-11  on  short,  finely-spaced  fins,  but  the 
accuracy  dropped  sharply  for  larger  fins  spaced  further 
apart.  This  was  apparently  due  to  the  assumed  linear  pres¬ 
sure  gradient  on  the  fin  surface,  and  breakdown  of  the  model 
as  gravity  forces  became  dominant.  The  use  of  equation 
(2.28)  should  be  restricted  to  fin  densities  from  1200  to 
IhOO  fins/m,  and  fin  heights  of  less  than  1  mm. 

In  1984,  Honda  et  al.  [24]  developed  a  model  to  predict 
the  heat-transfer  coefficient  that  took  surface-tension 
effects  and  the  non- isothermal  behavior  of  fins  into 
account.  In  order  to  obtain  a  solution,  however,  itera-ive 
finite  difference  techniques  that  required  a  considerable 
amount  of  computer  time  were  used.  Nonetheless,  they  showed 
that  this  model  was  able  to  predict  the  average  heat- 
transfer  coefficient  for  most  available  experimental  data 
(which  included  11  fluids,  and  22  tubes)  within  20  percent. 
Predictions  involving  the  condensation  of  steam  provided  the 
largest  errors  (as  much  as  40^)  ,  due  in  part  to  the  high 
surface  tension  of  water.  This  reflects  the  fact  that 
surface-tension  effects  are  still  not  completely  accounted 
for.  Compared  to  the  other  models  discussed  thus  far, 
however,  the  model  of  Honda  et  al.  shows  considerable 
promise  in  spite  of  its  complexity. 

In  1983,  Adamek  [27]  defined  a  family  of  condensate 
surface  profiles  whose  curvatures  are  given  by: 
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(2.29) 


1  =  Sm  fC  +  1) 

vhere 

=  rotation  angle  of  normal  to  condensate  film  surface, 

Sj^  =  length  of  convex  surface  of  fin  condensate  film, 

s  =  coordinate  along  condensate  film  surface  from 
fin  tip. 


The  parameter  i  in  equation  (2.29)  characterizes  the  aspect 
ratio  of  the  fin  cross  section.  k  number  of  C  values  and 
their  corresponding  condensate  surface  profiles  are  shown  in 
Figure  2.3.  !^ith  the  above  information,  idamek  obtained  an 
expression  for  the  heat-transfer  coefficient  of  the  fin 
surface  as  shown  below: 


=  2.140 
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(2.  30) 


These  fin  profiles  were  used  by  Hebb  et  al.  [25]  and  Rudy 
et  al.  [26]  in  1985  to  develop  a  model  for  condensation  on 
horizontal  finned  tubes.  This  model  divided  the  tube  into 
flooded,  and  unflooded  regions,  and  further  divided  the 
unflooded  region  into  tube  (area  between  fins)  and  fin 
areas.  To  predict  heat  transfer  from  the  unflooded  fin 
surfaces,  they  used  the  heat-transfer  coefficient  proposed 
by  Adamek.  The  Adamek  fin  profile  used  to  approximate  their 
trapezoidal  fins  is  shown  in  Figure  2.4.  The  Nusselt 


Figure  2.3  idaaek  [27]  Condensate  Surface  Profiles 


Figure  2.4  Fin  Geosetry  for  the  Rebb  et  al. 


Hodel  [22] 


was 


equation,  written  in  terms  of  the  film  Eeynolds  number, 
used  to  predict  heat  transfer  from  the  unflooded  tube 
surface  between  fins.  Writing  the  Nusselt  equation  this  way 
takes  into  account  the  additional  condensate  thickness  that 
results  from  the  drainage  of  condensate  formed  on  the  fin 
surface.  For  the  flooded  region,  a  two-dimensional  computer 
code  was  used  to  solve  for  the  heat  flux  into  the  tube-side 
coolant  ( »  which  was  compared  to  the  limiting  case  of 
heat  flux  (q^]^)  assuming  zero  fin  thickness.  This  ratio 
(41  =  qb2/'5bl  ^  used  along  with  a  linear  temperature 

profile  across  the  condensate  film  to  establish  a  heat- 
transfer  coefficient  for  this  region.  The  resulting  average 
heat- transfer  coefficient  is  given  below: 


h  =  h  n  = 
ow  o 


=  'h  ^  +  n  h.  1  +  h.  Ji; 

”  Af-f  ^  Apf  V  ^  IT 

2 

V  =  _ /-L..,.-..  Pe, 


h  =  4  IL 
^  e 


n  =  1  -  n  -  n't 


(2.31) 


(2.  32) 


(2.33) 


(2.  34) 


where  h ^  is  calculated  using  equation  (2.30),  and 


A gj  =  smooth  tube  area  based  on  the  fin  diameter. 


Re,  =  film  Reynolds  number. 


41 


'p 


ratio  of  tube  side  beat  flux  with  fins  to  tube 
side  heat  flux  with  fins  of  zero  thickness. 


A  letailed  step-by-step  procedure,  recommended  for  obtaining 
a  solution  of  equation  (2.31),  is  provided  in  Appendix  A. 

To  allow  for  heat  transfer  from  the  fin  tip,  and  to 
account  for  decreased  heat  transfer  due  to  the  thick  conden¬ 
sate  film  in  the  trough  area,  Webb  et  al.  replaced  ^  with 
Sj„  +  t/2  -  6^,  where  6^  is  the  average  condensate  film  thick¬ 
ness  in  the  trough  area  in  the  unflooded  region: 


(2.  35)  . 


Since  the  effects  due  to  t/2  and  6^  in  the  modified  tend 
to  cancel  each  other,  Webb  et  al,  reported  only  a  negligible 
change  in  S^,  This  model  was  shown  to  predict  the  heat- 
transfer  coefficient  for  R-11  condensing  on  horizontal 
finned  tubes  with  fin  densities  of  748,  1024,  and  1378 
fins/m  within  20  percent. 


Ill-  DESCRIPTION  OF  TEST  APPARATDS 

A.  TEST  APPARATUS 

The  same  test  apparatus  used  by  Seorgiadis  [7]  was  used 
for  this  investigation.  A  schematic  of  this  apparatus  is 
shown  in  Figure  3.1.  Steam  was  generated  in  the  boiler, 
which  consisted  of  a  304.8  mm  (12  in.)  Pyrex  glass  section 
containing  ten  4000-Ratt,  480-Volt  Fatlow  immersion  heaters. 
After  passing  through  a  304.8  mm  (12  in.)  to  152.4  mm 
(6  in.)  reducing  section,  the  steam  flowed  upward  through  a 
2.44  m  (8  ft.)  long  section  of  Pyrex  glass  piping.  At  this 
point,  a  180-degree  bend  in  the  piping  re-directed  the  steam 
downward,  where  after  1.52  m  (5  ft.)  of  Pyrex  glass  piping, 
it  entered  the  stainless  steel  test  section  illustrated  by 
Figure  3.2.  The  test  tube  was  mounted  horizontally  in  the 
test  section  behind  a  view  port,  which  allowed  visual  obser¬ 
vation  of  the  condensation  process.  Steam  not  condensed  by 
the  test  tube  continued  downward  to  the  auxiliary  condenser, 
where  it  was  condensed  by  two  9.5  mm  (3/8  in.)  diameter 
water-cooled  copper  tubes  helically  coiled  to  a  height  of 
457  mm  (18  in.).  The  condensate  was  then  returned  to  the 
boiler  by  gravity. 

Filtered  tap  water  was  used  to  cool  the  test  tube  on  a 
once-through  basis.  This  water  was  first  collected  in  a 
large  sump  with  a  capacity  of  about  0.4  cubic  meters  (Figure 
3.3),  then  pumped  through  a  flow  meter,  and  the  test  tube  by 
two  centrifugal  pumps  connected  in  series.  A  valve  on  the 
discharge  side  of  the  second  pump  allowed  the  velocity  of 
water  flowing  through  the  test  tube  to  be  varied  from  0  to 
4.4  m/s  (14.4  ft/sec).  The  auxiliary  condenser  was  cooled 
by  a  continuous  supply  of  tap  water,  which  was  throttled  to 
control  the  internal  pressure  of  the  test  apparatus. 
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The  system  used  to  remove  Doa-condensing  gases  is  shown 
in  Figure  3.3.  The  vacuum  pump  continually  drew  a  sample 
mixture  of  any  possible  non-condensing  gases  and  moisture 
from  the  auxiliary  condenser.  This  mixture  was  then  passed 
through  a  heat  exchanger,  where  the  moisture  was  condensed 
before  it  entered  the  sump,  using  the  same  filtered  tap 
water  described  earlier. 

B.  INSTBUHENTATIOM 

The  electrical  power  input  to  the  boiler  immersion 
heaters  was  controlled  by  a  panel-mounted  potentiometer.  A 
converter  with  an  input  voltage  of  440  VAC  generated  a 
signal  which  was  fed  to  the  data  acquisition  system  to 
calculate  the  boiler  input  power.  A  more  detailed  descrip¬ 
tion  of  the  boiler  power  supply  is  provided  by  Poole  [6]. 

The  internal  system  pressure  was  measured  manually  using 
a  0-tube,  mercury- in-glass  manometer  graduated  in  millime¬ 
ters.  Steam,  condensate  return,  and  ambient  temperatures 
were  measured  using  calibrated  copper-constantan  thermocou¬ 
ples.  These  thermocouples  had  an  accuracy  within  ±  0.1  K. 
The  cooling  water  temperature  rise  was  measured  by  two 
Hewlett-Packard  (HP)  2804A  quartz  thermometers,  along  with  a 
10- junction,  series-connected  copper-constantan  thermopile 
as  a  backup.  Throughout  this  investigation,  the  quartz 
thermometers  and  the  thermopile  agreed  to  within  ±  0.03  K. 

C.  DATA  ACQUISITIOH  STSTEH 

To  monitor  the  system  temperatures  and  boiler  input 
power  (using  the  converter  signal) ,  an  HP  9826A  computer  was 
used  to  control  an  HP  3497A  Data  Acquisition  System.  Eaw 
data  were  processed  immediately  and,  at  the  same  time, 
stored  on  diskett  for  reprocessing  at  a  later  time. 
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D.  TUBES  TESTED 


For  this  thesis  effort,  eighteen  copper  tubes,  and  a 
stairless  steel  tube  were  manufactured.  They  are  presented 
in  Table  I.  The  first  eight  tubes  with  rectangularly-shaped 
fins  were  tested  in  order  to  complete  the  seguence  of  data 
initiated  by  Seorgiadis  [7].  This  completed  the  systematic 
variation  of  fin  spacing,  thickness,  and  height  so  that  the 
optimum  dimensions  could  be  determined.  Figure  3.4  shows  a 
photograph  of  the  four  tubes  with  fin  thickness  of  1.0  mm, 
fin  height  of  0.5  mm,  and  spacings  of  1.0,  1.5,  2.0,  and  4.0 
mm  (tube  numbers  26  to  29).  Also  shown  is  the  insert  used 
to  enhance  the  inside  heat- transfer  coefficient,  which  will 
be  discussed  in  more  detail  in  the  next  chapter.  Figure 
3.5(a)  shows  a  cross  section  of  a  rectangular  fin  with  a  fin 
height  of  1-.  5  mm,  a  fin  spacing  of  1.5  mm,  and  a  fin  thick¬ 
ness  of  1.0  mm.  In  addition  to  the  tubes  with 
rectangularly-shaped  fins,  six  tubes  were  tested  to  study 
the  effect  of  fin  shape  on  the  heat-transfer  performance. 
Tubes  number  40  and  42  were  made  of  stainless  steel  to  test 
the  effect  of  fin-metal  thermal  conductivity  on  the  heat- 
transfer  performance.  Due  to  limited  machine-shop  facili¬ 
ties  available,  the  ’’parabolic"  fins  on  two  of  the  tubes 
were  not  truly  parabolic,  as  shown  in  Figure  3.5(b).  A 
cross  section  of  a  commercially-available  finned  tube 
produced  by  the  Wolverine  Division  of  Onited  Oil  Products  is 
also  shown  in  Figure  3.5(c)  for  comparison.  The  smooth 
tubes  were  tested  in  order  to  determine  the  water-side  heat- 
transfer  coefficients,  and  to  serve  as  a  comparison  for  the 
finned  tubes. 

E.  VACOOH  IITEGHITI 

For  any  condensation  experiment,  especially  when  oper¬ 
ating  under  vacuum  conditions,  a  leak- free  apparatus  is  of 
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vital  importance.  The  presence  of  any  leaks  would  allow  the 
admission  of  non-condensing  gases,  whose  deleterious  effects 
would  invalidate  any  data  being  collected.  To  prevent  this 
problem,  vacuum  integrity  tests  were  performed  on  the  test 
apparatus  both  before  and  after  the  data  acquisition  phase 
of  this  thesis  effort  (a  time  period  of  about  four  months) . 
The  first  test  revealed  a  leak  rate  represented  by  a  pres¬ 
sure  rise  of  1.0  mmHg  over  a  period  of  five  days,  while  the 
pressure  rise  during  the  final  test  was  1.0  mmHg  over  a 
24-hour  period.  These  tests  were  performed  at  an  absolute 
pressure  of  about  85  mmHg.  This  negligible  leak  rate, 
together  with  continuous  venting  of  the  test  apparatus 
during  data  runs,  eliminated  the  deleterious  effects  of 
non-condensing  gases  on  the  heat-transfer  results  reported 
in  this  thesis. 


IT.  DATA  CQLLECTIOH  ABD  BEDgCTIOH 
A.  SYSTEM  OPERATION 

In  order  to  ensure  filmwise  condensation,  all  tubes  were 
treated  with  a  solution  containing  equal  parts  of  sodium 
hydroxide  and  ethyl  alcohol.  The  tubes  were  placed  in  a 
steam  bath  and  several  layers  of  the  solution  were  applied 
with  a  tooth  brush  every  10  minutes.  After  about  one  hour, 
the  tube  surfaces  were  completely  blackened.  The  same 
procedure  was  repeated  prior  to  each  data  run,  but  the  tube 
was  only  steamed  for  tec  to  fifteen  minutes.  Following  the 
steam  bath,  the  tubes  were  rinsed  with  tap  water  before 
Installing  them  in  the  test  apparatus.  Georgiadis  [7] 
discusses  this  procedure  in  more  detail. 

The  test  apparatus  was  brought  to  operating  pressure  and 
temperature  using  the  procedures  provided  by  Georgiadis  [7]. 
Steady-state  conditions  were  assumed  to  occur  once  the  test 
section  steam  temperatures,  as  well  as  the  cooling  water 
inlet  and  outlet  temperatures  stabilized.  Two  sets  of  data 
were  then  taken  for  each  of  the  heat  fluxes,  which  were 
determined  by  the  flow  rate  of  cooling  water  through  the 
test  tube.  Starting  at  a  flow  rate  of  80  percent  (4.4U  ra/s 
for  19  mm  O.D.  tubes,  and  3.92  m/s  for  13.7  and  15.8  mm  C.D. 
tubes),  flow  rates  were  reduced  to  70,  60,  45,  35,  26,  and 
20  percent,  then  brought  back  up  to  55  and  90  percent. 
These  settings  were  selected  since  they  provided  nearly 
equally-spaced  heat  flux  values.  It  should  be  noted  that 
continuous  adjustments  of  cooling  water  flow  through  the 
auxiliary  condenser  were  required  to  maintain  system  pres¬ 
sure  as  flow  rates  through  the  test  tube  were  changed. 


visual  observations  of  condensation  taking  place  on  the 
test  tube  were  made  on  a  regular  basis  to  ensure  uniform 
filmwise  condensation.  A  more  reliable  indication  of  this, 
however,  was  the  trend  of  data  taken  as  the  cooling  water 
flow  rates  were  increased  back  up  to  80  percent.  If  the 
final  value  of  the  steam-side  heat- transfer  coefficient  at 
80  percent  was  significantly  higher  than  those  at  the  begin¬ 
ning  of  the  data  run,  dropwise  condensation  was  assumed  to 
have  occurred,  and  the  data  were  disregarded.  All  data 
presented  in  this  thesis  displayed  less  than  555  disagreement 
in  the  steam-side  heat-transfer  coefficient  between  initial 
and  final  data  sets. 

B.  THE  DBOPHISE  CONDENSATION  PBOBLEIi! 

The  purpose  of  this  thesis  was  to  take  strictly  filmwise 
condensation  data,  using  primarily  copper  tubes.  Due  to  the 
poor  wettability  of  copper  with  water,  especially  when  even 
minute  amounts  of  contamination  are  present,  there  is  a 
tendency  to  condense  steam  in  the  partial  dropwise  mode. 
Dropwise  condensation  is  far  more  effective  than  filmwise 
condensation;  so,  its  presence  can  lead  to  large  errors  in 
the  data. 

By  following  procedures  set  forth  by  Georgiadis  [7],  the 
dropwise  problem  was  minimized.  During  the  initial  testing 
phase  of  this  thesis,  however,  a  small  amount  of  dropwise 
condensation  was  encountered.  This  problem  was  soon  traced 
to  the  boiler  water,  which  had  not  been  changed  for  some 
time.  By  changing  the  water  after  every  three  or  four  runs, 
the  problem  was  eliminated. 

C.  STEAH  VELOCITY  LIHITATIONS 

A  major  assumption  in  the  Nusselt  theory  is  that  vapor 
velocity  does  not  induce  shear  forces  on  the  vapor-liguid 


interface.  Since  system  pressure  in  the  test  apparatus  was 
controlled  by  condensing  steam  in  the  auxiliary  condenser,  a 
reasonable  steam  velocity  was  always  present.  Since  steam 
velocity  (i.e. ,  steam  flow  rate)  depends  on  boiler  power 
input,  decreasing  the  power  results  in  lover  vapor  veloci¬ 
ties.  Dnder  vacuum  conditions,  (about  85  mmHg)  a  maximum 
velocity  of  about  8  m/s  was  possible  at  full  boiler  power 
(36  kF) .  However,  at  these  high  velocities,  the  entire 
apparatus  experienced  considerable  vibration.  As  discussed 
by  Georgiadis  [7],  a  steam  velocity  of  2  m/s  under  vacuum 
provided  the  most  stable  operating  conditions.  This  corre¬ 
sponding  boiler  power  (about  9  kH)  was  low  enough  to  mini¬ 
mize  system  vibrations,  yet  high  enough  to  avoid  rapid 
fluctuations  in  condenser  pressure  due  to  the  intermittent 
break-up  of  vapor  bubbles  in  the  boiler.  To  allow  a  suffi¬ 
cient  amount  of  steam  to  reach  the  auxiliary  condenser  when 
operating  at  atmospheric  pressure,  the  maximum  boiler  power 
was  required.  This  resulted  in  a  steam  velocity  of  approxi¬ 
mately  1  m/s. 

D.  DATA  REDDCTIOH 

All  the  programs,  property  functions,  and  calibration 
data  used  for  this  thesis  effort  were  essentialy  the  same  as 
those  used  by  Georgiadis  [7]  and  Poole  [6].  Because  six  of 
the  tubes  (tube  numbers  39  thru  44)  tested  during  this 
effort  had  inside  diameters  different  than  any  tube  tested 
previously,  and  two  (tubes  number  40,  and  42)  were  manufac¬ 
tured  from  stainless  steel,  some  minor  modifications  were 
necessary  for  the  data  reduction  program  (program  DRP6)  . 
These  modifications  included  options  for  different  tube 
diameters,  tube  inserts,  thermal  conductivity,  mixing- 
chamber  calibrations,  and  leading  constants  for  the  inside 
heat-transfer  coefficients. 


The  mixing- chamber  (see  Figure  3.2)  calibration  was 
required  to  account  for  the  temperature  rise  that  occurred 
across  the  mixing  chamber,  due  primarily  to  turbulent 
viscous  dissipation.  The  purpose  of  the  mixing  chamber  is 
to  ensure  a  uniform  temperature  at  the  cooling  water 
discharge,  so  that  accurate  temperature  measurements  can  be 
made.  Since  tubes  with  different  inside  diameters, 
requiring  different  inserts,  result  in  different  pressure 
drops,  calibrations  had  to  be  performed  for  each  inside 
diameter. 

A  modified  Wilson  Plot  program  [7]  was  used  to  compute 
the  leading  constants  for  the  Sieder-Tate  equation  used  to 
determine  the  water-side  heat-transfer  coefficient.  Data 
taken  on  smooth  tubes  were  used  to  determine  the  leading 
constants,  which  were  then  used  to  reduce  data  taken  on 
similar  tubes  with  identical  internal  dimensions. 


7.  BESOLTS  MS  DISC  OSS  ION 

A.  INTRODOCTION 

During  this  thesis  effort,  a  number  of  data  runs  were 
made  using  the  procedures  described  in  Chapter  IV.  For  both 
vacuum  and  atmospheric  conditions,  every  tube  was  tested  at 
least  three  times,  on  different  days,  to  ensure  the  data 
were  repeatable.  Complete  filmwise  condensation  was  main¬ 
tained,  and  the  computed  non-condensing  gas  concentration 
was  held  below  1.0  percent  for  all  accepted  data  runs.  It 
was  estimated  that  an  error  of  ±  0.2  K  in  the  steam  tempera¬ 
ture,  or  an  error  of  ±  1.0  mmEg  in  the  system  pressure  would 
lead  to  an  error  of  ±  1.2  percent  in  the  non-condensing  gas 
concentration.  This  shows  that  the  computed  non-condensing 
gas  concentration  was  zero  to  within  the  accuracy  of  temper¬ 
ature  and  pressure  measurements.  As  discussed  earlier  in 
Chapter  III,  the  test  apparatus  would  allow  only  a  negli¬ 
gible  amount  of  non-condensing  gas  to  be  present  in  the 
apparatus.  The  non-condensing  gas  concentration  was 
computed  for  the  system  temperature  and  pressure  during 
every  data  run  to  ensure  no  major  leaks  developed.  A 
summary  of  tubes  tested  between  Georgiadis  [7]  and  this 
thesis  effort  and  the  resulting  enhancements  are  provided  in 
Table  II. 

B.  NATEB-SIDE  HEAT-TBANSFBB  COEFFICIENTS 

The  water-side  heat-transfer  coefficients,  enhanced  by 
the  use  of  an  insert,  were  obtained  by  Georgiadis  [7]  for 
the  tubes  with  an  inside  diameter  of  12.7  mm.  Both  a 
’’direct  method,”  and  a  ’’modified  Wilson  method”  were  used, 
and  are  discussed  below. 


TABf.K  II 


1.  The  Direct  Method 


For  the  "direct  method,"  a  smooth  tube  was  manufac¬ 
tured  with  six  wall  thermocouples  inssrted  around  its  perim¬ 
eter.  During  data  runs  under  vacuum  and  atmospheric 
conditions,  readings  from  these  six  thermocouples  were 
obtained  and  averaged,  in  order  to  determine  an  average  tube 
wall  temperature  for  each  heat  flux.  It  should  be  noted 
that  the  wall  temperature  was  guite  non-uniform  around  the 
tube  perimeter.  Georgiadis  [7]  reported  up  to  an  18-K  drop 
in  the  wall  temperature  from  the  top  to  the  bottom  of  the 
tube.  Evaluation  of  the  uncertainty  resulting  from  this 
temperature  variation  was  considered  beyond  the  scope  of  the 
investigation.  The  average  wall  temperature  was  then  used 
along  with  the  cooling  water  inlet  and  outlet  temperatures 
to  determine  the  log-mean-temperature  difference  (LMTD)  . 
With  values  for  heat  flux  and  LMTD  known,  the  value  of  the 
water-side  heat- transfer  coefficient  was  obtained.  This,  in 
turn,  was  used  to  determine  the  leading  coefficient  for  the 
Sieder-Tate-type  equation  shown  below: 


Nn  =  a  Re°'®  T>r 


1/3  fn 
_l‘w_ 


(5.1) 


A  more  detailed  step-wise  solution  procedure  is  given  by 
Georgiadis  [7].  For  the  previously-mentioned  tubes  of 
15.8  mm  and  19.05  mm  O.D.,  Georgiadis  determined  the  coeffi¬ 
cient  C  to  be  0.0615.  The  constant  B  =  26.4  is  an  addi¬ 
tional  parameter  found  to  improve  the  fit  of  equation  (5.1). 
Dote  that  C  =  0.0635  is  about  2.5  times  greater  than  the 
well-known  Sieder-Tate  coefficient  of  0.027  [28]  for  long 
tubes  with  smooth  inner  walls.  This  is  due  primarily  to  the 
insert  used  to  enhance  the  water-side  heat-transfer 
coefficient. 


A  complete  description  of  this  method  is  provided  in 
Appendix  B.  Briefly,  this  method  uses  a  Sieder-Tate-type 
equation  for  the  water-side  heat-transfer  coefficient,  and 
an  equation  (discussed  below)  for  the  steam-side  heat- 
transfer  coefficient.  Both  equations  contain  leading  coef¬ 
ficients  that  must  be  determined  iteratively  while  being 
used  in  a  linear  equation  to  generate  the  Wilson  plot.  The 
iterations  were  performed  on  the  HP  9826A  computer  mentioned 
previously  in  Chapter  III,  using  program  •* WILSON 6. '• 

The  "modified  Wilson  method"  used  during  this  inves¬ 
tigation  consisted  of  one  minor  modification  from  that  of 
Georgiadis.  For  this  thesis,  a  Fujii-type  [29]  equation  was 
used  instead  af  the  Nusselt-type  equation  for  the  steam-side 
coefficient.  This  modification  was  made  to  account  for  the 
small  steam  velocity  (about  1  m/s)  that  was  present  during 
the  runs.  For  this  reason,  use  of  the  Fujii-type  equation 
appears  to  be  more  accurate  than  the  use  of  the  Nusselt-type 
equation.  This  modification  resulted  in  a  slightly  higher 
value  (up  to  3  percent)  for  the  leading  coefficient  C.  The 
program  "FI1S0N6"  described  earlier,  allows  an  option  for 
selecting  either  the  Fujii-type  or  Nusselt-type  equation  for 
the  steam-side  coefficient. 

Using  a  Nusselt-type  equation  for  the  steam-side 
coefficient,  Georgiadis  [7]  found  the  leading  coefficient 
(for  equation  (5.1))  C  to  be  0.071,  with  the  B  value  set 
equal  to  zero.  This  C  value  is  about  10  percent  higher  than 
the  value  (C  =  0.0635)  obtained  by  the  "direct  method,"  but 
since  the  "direct  method"  is  generally  felt  to  be  more  reli¬ 
able,  the  values  of  0.0635  and  26.4  were  used  for  the 
constants  C  and  B  respectively  for  this  thesis  effort  with 
tubes  of  12.7  mm  inside  diameter  (see  Table  I). 


3.  Water-Side  Coefficients  For  Thin- Wall  Tubes 

Due  to  the  thin  (i.e.,  0.5  mm  thickness)  tube  vails 
for  tubes  39  thru  42,  it  was  not  possible  to  manufacture  an 
Instrumented  tube.  For  this  reason,  the  "modified  Wilson 
method,"  together  with  the  Fujii-type  equation  for  the 
steam-side  coefficient,  was  used  to  determine  the  leading 
coefficient  for  equation  (5.1)  for  both  the  copper  tube 
(tube  number  41)  ,  and  the  stainless  steel  tube  (tube  number 
42)  .  For  the  copper  tube,  a  value  of  C  =  0.0756  was 
obtained,  while  for  the  stainless  steel  tube,  a  value  of 
C  =  0.0688  was  obtained.  Although  the  inner  tube  diameters 
for  both  tubes  are  the  same  (i.e.,  13.5  mm),  the  much  higher 
thermal  conductivity  of  copper,  which  results  in  circumfer¬ 
ential  heat  conduction  through  the  tube  wall,  leads  to  its 
higher  value  for  C.  A  summary  of  the  leading  coefficients 
used  in  equation  (5.1)  for  all  the  tubes  is  presented  in 
Table  III. 

C.  EEPEATABILITI  OP  DATA 

Since  the  reliability  of  data  taken  was  of  vital  impor¬ 
tance,  all  data  runs  were  repeated,  as  Georgiadis  [7]  did, 
at  least  three  times  on  different  days.  The  computed  steam- 
side  coefficients  for  similar  conditions  (i.e.,  same  tube 
and  the  same  operating  conditions)  from  different  days 
agreed  to  within  ±  5  percent.  Additionally,  data  runs  were 
performed  on  two  finned  tubes  (tubes  number  6  and  17),  under 
similar  conditions  to  verify  the  repeatability  with  data 
taken  by  Georgiadis.  Figure  5.1  shows  that  the  experimental 
steam-side  heat-transfer  coefficients  of  Georgiadis  and 
those  obtained  during  this  investigation  agree  to  within  ±  5 
percent.  The  curves  shown  in  this  figure  (and  all  other 
figures  of  a  similar  nature  that  follow)  are  the  least- 
sguares-fit  curves  according  to  the  following  equation: 


CoBparisoo  of  Finoed  Tube  Data  with 


where  a  and  b  are  experimentally  determined  constants. 
These  curves  were  selected  instead  of  those  based  on  the 
steam-side  coefficient  versus  heat  flux  since  the  uncer¬ 
tainty  band  on  AT  is  more  uniform  than  the  uncertainty  band 
based  on  the  steam-side  coefficient.  A  brief  discussion  of 
the  uncertainty  bands  is  presented  in  Appendix  D,  along  with 
a  listing  of  the  uncertainty  analysis  program  ”0NA6''  and  a 
few  sample  runs. 


D.  EFFECTS  OP  FIH  SPACING  AND  FIN  HEIGHT  ON  PEfiFOBHAHCE 

Data  were  taken  on  eight  tubes  (see  Table  II)  with  fins 
of  rectangular  shape  (tubes  number  26  thru  33)  under  vacuum 
(approximately  85  mmllg)  and  atmospheric  conditions.  Four 
tubes  had  a  fin  height  of  0.5  mm,  and  fin  spacings  of 
1.0,  1.5,  2.0,  and  h.O  mm,  while  the  remaining  four  tubes 
had  a  fin  height  of  1.5  mm,  and  the  same  seguence  of  fin 
spacings.  Georgiadis  [7]  reported  data  on  rectangularly- 
shaped  fins  with  the  same  fin  spacings,  but  with  fin  heights 
of  1.0  mm,  and  2.0  mm.  A  complete  listing  of  tubes  tested 
by  Georgiadis  is  also  presented  in  Table  II  (tubes  1  thru 
25)  . 

Data  reductions  were  performed  on  the  KP  9326A  computer, 
using  a  program  named  '’DEP6.  "  This  program  was  essentially 
the  same  one  used  by  Georgiadis  (''DRP4”)  ,  with  modifications 
to  allow  for  variations  in  tube  diameters,  and  tube 
material.  Data  obtained  for  all  tubes  tested  during  this 
thesis  effort  (tubes  number  26  thru  44)  are  presented  in 
Appendix  C. 


1. 


Effects  of  Fin  Spacing 

This  section  presents  results  showing  the  variations 
of  the  steam-side  heat-transfer  coefficient  with  heat  flux 
having  fin  spacing  as  a  parameter.  Figures  5.2  and  5.3 
present  data  for  tubes  with  a  fin  height  of  0.5  mm  under 
vacuum  and  at  atmospheric  pressure  respectively.  The  smooth 
tube  data  and  the  theoretical  prediction  of  Nusselt  are  are 
included  for  comparison.  In  a  similar  way.  Figures  5.4  and 

5.5  present  data  for  tubes  with  a  fin  height  of  1.5  mm. 
These  figures  show  that  the  best  heat-transfer  performance 
occurs  for  the  tubes  with  a  fin  spacing  of  2.0  mm,  while  the 
worst  performance  occurs  for  the  tubes  with  a  fin  spacing  of 
4.0  mm.  Note  that  these  trends  are  the  same  both  under 
vacuum  and  at  atmospheric  conditions. 

The  comparison  of  finned  tubes  is  made  through  the 
enhancement  ratio.  This  ratio  is  defined  as  the  steam-side 
heat-transfer  coefficient  of  a  finned  tube  to  that  of  the 
smooth  tube  (same  diameter  as  the  finned  tube  root  diameter) 
at  the  same  heat  flux.  As  was  shown  by  Georgiadis  [7],  the 
uncertainty  in  calculating  the  steam-side  heat- transfer 
coefficient  increases  with  decreasing  heat  flux;  so,  heat 
flux  values  of  0.35  and  1.0  MW/m^  were  chosen  for  vacuum  and 
atmospheric  conditions  respectively.  Cross  plots  of  the 
enhancement  ratio  versus  fin  spacing  are  shown  in  Figures 

5.6  and  5.7  for  fin  heights  of  0.5,  1.0,  1.5.  and  2.0  mm 
(curves  representing  e  =  1.0  and  2.0  mm  are  from  Georgiadis 
[7]).  Enhancements  as  high  as  4.8  and  6.4  can  be  seen  for 
vacuum  and  atmospheric  conditions  respectively.  Once  again, 
fin  spacings  of  approximately  2.0  mm  for  the  tubes  with  fin 
heights  of  0.5  and  1.5  mm  are  shown  to  provide  the  optimum 
performance.  Data  provided  by  Georgiadis  ^7]  for  fin 
heights  of  1.0  and  2.0  mm,  however,  demonstrate  an  optimum 
heat-transfer  performance  for  fin  spacings  of  1.5  mm.  As  a 
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5.7  Cross  Plot  Showing  Best  Pin  Spacing  Aaong  Tabes 
ith  e  0.5,  1.0,  1.5,  and  2.0  nm  (Xta.  kuqs)  . 


result,  there  3oes  not  appear  to  be  a  clear  trend  for  the 
optimum  fin  spacing  as  a  function  of  fin  height. 

On  one  hand,  the  difference  in  performance  between 
the  s  1. 5  and  2.0  mm  finned  tubes,  as  well  as  the  s  -  0.5 
and  1.0  mm  finned  tubes  is  less  than  the  experimental  uncer¬ 
tainty  (±  15  percent),  so  that  the  optimum  spacing  may  be 
somewhere  between  1.5  and  2.0  mm.  On  the  other  hand,  each 
data  point  in  these  figures  has  been  repeated  within  5 
percent  (i.e. ,  computed  steam-side  heat-transfer  coeffi¬ 
cient)  on  different  days.  Also,  as  discussed  earlier,  the 
data  taken  on  tubes  6  and  17  by  Georgiadis  were  repeated 
within  3  percent  during  this  investigation.  Therefore,  the 
slightly  disagreeing  trends  shown  by  Georgiadis  and  by  this 
investigator  (see  Figures  5.6  and  5.7)  may  be  justified  by 
the  very  complex  nature  of  the  condensation  phenomenon  on 
finned  tubes. 

The  presence  of  an  optimum  fin  spacing,  however,  is 
easy  to  understand.  The  1,0-mm  fin  spacing,  which  provides 
the  largest  increase  in  surface  area,  performed  worse  than 
the  2.0-mm  fin  spacing  due  to  the  relatively  large  conden¬ 
sate  retention  angle  with  the  1.0  mm  fin  spacing.  Visual 
observations  showed  that  the  condensate  retention  angle  was 
around  110  degrees  for  the  tube  with  s  =  1.0  mm,  while  it 
was  about  40  degrees  for  the  tube  with  s  =  2.0  mm.  The 
additional  thermal  resistance  induced  by  this  layer  of 
condensate  overpowers  the  benefit  gained  from  increased 
surface  area,  so  the  heat-transf er  performance  is  reduced. 
As  fin  spacing  increases,  the  condensate  retention  angle 
decreases,  as  does  the  thickness  of  the  condensate  film 
between  fins  on  the  upper  portion  of  the  tube.  Beyond  fin 
spacings  of  2.0  mm,  the  decrease  in  surface  area  tends  to 
overshadow  the  other  effects,  resulting  in  decreased  heat- 
transfer  performance  compared  to  the  tube  with  the  "optimum" 
spacing. 


2*  Effects  of  Fin  Height 

As  mentioned  earlier,  fin  heights  of  0.5  and  1.5  mm 
were  examined  during  this  thesis  effort,  and  data  for  the 
tubes  with  fin  heights  of  1.0  and  2.0  mm  were  available  from 
Georgiadis  [7].  Table  II  shows  that  the  general  trend  is  an 
increase  of  enhancement  ratio  as  fin  height  increases,  which 
results  from  the  increase  of  surface  area  provided  by  the 
higher  fins. 

To  study  the  enhancement  obtained  beyond  the 
increase  in  surface  area  due  to  finning,  two  additional 
columns  are  provided  in  Table  II.  These  columns  represent 
enhancement  catio/area  ratio  (£o/Ar)  at  both  vacuum  and 
atmospheric  conditions.  As  can  be  seen,  Eo/Ar  values  range 
from  0.86  to  2.2  3  under  vacuum  conditions,  and  from  1.43  to 
3.20  at  atmospheric  pressure.  Figures  5.8  and  5.9  present 
plots  of  of  Eo/Ar  as  a  function  of  fin  spacing  under  vacuum 
conditions,  and  at  atmospheric  pressure  respectively.  The 
relatively  small  values  of  Eo/Ar  =  0.86  (vacuum)  and 
Eo/Ar  =  1.43  (atmospheric)  for  tube  number  15  are  due  to  the 
small  fin  spacing  which  results  in  a  completely  flooded 
tube.  On  the  other  hand,  the  values  of  Eo/Ar  =  2.23 
(vacuum)  for  tube  number  14,  and  Eo/Ar  =  3.20  (atmospheric) 
for  tube  number  28,  show  that  significant  enhancements 
beyond  the  area  ratio  are  obtainable  when  the  amount  of 
condensate  flooding  is  kept  small.  It  appears  that  surface- 
tension  forces  are  mostly  responsible  for  this  enhancement 
beyond  the  area  ratio. 

As  discussed  in  Chapter  II,  the  condensate  thinning 
provided  by  surface-tension  effects  is  related  to  the  pres¬ 
sure  gradient  from  the  fin  tip  to  the  fin  root.  For  short 
fins,  this  pressure  gradient  is  relatively  larger  than  for 
high  fins,  resulting  in  improved  thinning  of  the  condensate 
film.  On  the  other  hand,  as  fin  height  increases,  surface 


Cross  Plot  Showing  Best  Eo/Ar  Along  Tabes 
0.5r  1*0,  1>5,  and  2.0  ■■  (Ata.  Bans). 


area  increases  and  the  condensate  retention  angle  decreases. 
Both  of  these  effects  lead  to  increased  heat- transfer 
perf oricance,  but  when  combined  with  the  poorer  thinning 
effect,  the  Eo/Ar  ratio  decreases  somewhat  for  the  larger 
fin  heights. 

E.  EFFECT  OF  FIH  GECHETRT  ON  PEBFOEMANCE 

In  order  to  study  the  effect  of  fin  geometry,  data  were 
taken  on  tubes  with  fins  of  triangular,  trapezoidal,  and 
parabolic  shapes  (tubes  number  34,  35,  37,  38)  ,  as  well  as 
on  a  spirally- finned  tube  with  threads  of  triangular  shape 
(tube  number  36).  Fin  dimensions  (i.e.,  fin  spacing, 
height,  and  thickness)  for  these  tubes  were  chosen  so  that 
the  resulting  total  surface  areas  would  be  approximately 
equal  to  that  of  the  "optimum*'  tube  with  rectangularly- 
shaped  fins  and  a  fin  height  of  1.0  mm  (tube  number  6) .  In 
this  manner,  the  effects  of  increased  surface  area  were 
eliminated,  allowing  a  comparison  between  tubes  where  fin 
geometry  was  the  only  variable.  A  commercially-available 
"Kblverine"  tube  (tube  number  43),  whose  fin  geometry  is 
shown  in  Figure  3.5(c),  was  also  tasted.  Dimensions  for  the 
above  fins  are  given  in  Table  II. 

1 .  Effect  of  Fin  Shape 

This  section  examines  the  effect  that  fin  shape  has 
on  the  heat-transfer  performance  of  finned  tubes.  The 
performance  of  tubes  number  34,  35,  36,  37,  and  43  under 
vacuum  conditions  is  shown  in  Figure  5.10,  while  Figure  5.11 
depicts  their  performance  at  atmospheric  pressure.  The 
dashed  curve  in  these  figures  represents  a  least-squares  fit 
for  data  taken  on  the  "optimum"  tube  with  rectangularly- 
shaped  fins  and  a  fin  height  of  1.0  mm  (tube  number  6).  In 
both  cases,  the  "Wolverine"  tube  (tube  number  43)  gave  the 


Figure  5.10  Tariatioo  of  Heat-Transfer  Coefficient  for  Tubes 
with  Fins  of  Different  Shapes  (Vac.  Runs). 


worst  performance.  This  was  due  primarily  to  the  small  fin 
spacing  (0.7  mm)  that  resulted  in  a  completely  flooded  tube. 
Based  on  the  results  of  Georgiadis  [7],  the  small  fin  thick¬ 
ness  (0.3  mm)  of  the  Wolverine  tube  may  also  be  responsible 
for  its  poor  performance.  Georgiadis  showed  that  for  a  fin 
spacing  of  1.0  mm  the  tubes  with  a  fin  thickness  of  0.5  mm 
(tube  number  16)  performed  16  percent  and  6  percent  poorer, 
under  vacuum  and  at  atmospheric  pressure  respectively,  than 
the  tube  with  a  fin  thickness  of  1.0  mm  (tube  number  5). 
Eased  on  this  observation,  the  small  fin  thickness  is  also 
responsible  for  the  poorer  performance  of  the  "Wolverine" 
tube  in  comparison,  for  example,  to  the  tube  with  fins  of 
parabolic  shape  (tube  number  37),  which  was  also  completely 
flooded  by  condensate.  Therefore,  the  poor  performance  of 
this  tube  is  probably  due  more  to  increased  condensate 
flooding  and  smaller  fin  thickness  than  to  fin  shape.  On 
the  other  hand,  the  best  performance  was  observed  for  the 
spirally-finned  tube  (tube  number  36)  under  vacuum  condi¬ 
tions,  and  the  tube  with  parabolic  fins,  where  fin  spacing 
was  0.7  mm  (tube  number  37)  ,  at  atmospheric  pressure.  The 
fact  that  tube  number  37  performed  well  at  atmospheric  pres¬ 
sure,  but  poorer  than  the  spirally-finned  tube  under  vacuum 
conditions  was  somewhat  unexpected.  It  must  be  kept  in  mind 
that  these  data  runs  were  repeated  at  least  three  times  on 
different  days,  thus  confirming  this  trend.  One  possible 
explanation  centers  on  variable  fluid  properties,  as 
discussed  below. 

Under  vacuum  conditions,  the  condensate  retention 
angle  is  significantly  smaller  for  tube  number  36  (about  30 
degrees)  than  for  tube  number  37  (completely  flooded) ,  and 
the  saturation  temperature  is  about  50  K  lower.  Therefore, 
the  relatively  high  heat-transfer  coefficient  of  the 
spirally-finned  tube  is  due  mainly  to  its  small  condensate 
retention  angle.  When  the  condensing  temperature  is  raised. 


both  surface  tension  and  viscosity  decrease.  For  example, 
the  surface  tension  decreases  by  about  10  percent,  and  the 
viscosity  decreases  by  about  45  percent  when  the  temperature 
is  redsed  from  40  oc  to  90  oc  (assumed  approximate  film 
temperatures).  Decreased  viscosity  should  result  in  higher 
heat-transf er  coefficients  at  atmospheric  pressure  than 
under  vacuum  conditions  for  a  given  tube.  Further,  it  may 
be  possible  that  the  effect  of  viscosity  is  dependent  on  the 
fin  shape.  Dnfortunately ,  no  eguations  exist  to  support 
this  hypothesis.  The  effect  of  surface  tension,  on  the 
other  hand,  is  fairly  well  understood;  the  condensate  reten¬ 
tion  angle  decreases  with  decreasing  surface  tension. 
Therefore,  the  reduction  in  retention  angle  (from  vacuum  to 
atmospheric  pressure)  for  the  tube  with  "parabolic"  fins  is 
greater  than  the  reduction  for  the  spirally-finned  tube, 
since  the  latter  was  already  small  under  vacuum  conditions. 
Thus,  the  reduced  retention  angle  on  the  tube  with  "para¬ 
bolic"  fins  provides  an  increased  heat-transf er  performance, 
surpassing  the  performance  of  the  spirally-finned  tube. 

The  performance  of  the  remaining  two  tubes  (tubes 
number  34  and  35)  is  the  same  under  vacuum  conditions,  but 
at  atmospheric  pressure,  the  tube  with  triangularly-shaped 
fins  (tube  number  34)  outperformed  the  tube  with  trapezoidal 
fins  (tube  number  3  5)  by  about  10  percent.  An  argument 
similar  to  the  one  above  may  explain  this  behavior. 


2.  Effects  on  Enhancement  Eatio 


Referring  to  Table  II,  it  can  be  seen  that  all  of 
the  above  tubes  provided  significant  heat-transfer  enhance¬ 


ments  over  smooth-tube  values. 


Mote  that  the  enhancement 


obtained  under  vacuum  conditions  for  tube  number  36  was  the 
same  as  obtained  for  the  "optimum,"  rectangularly-shaped 
finned  tube  with  the  same  fin  height  (tube  number  6) .  Under 
atmospheric  conditions,  however,  the  tube  with  fins  of 


parabolic  shape  (tube  number  37)  performed  significantly 
better  than  any  of  the  tubes  with  rectangular  fins  of  the 
same  height  (1.0  mm)  . 

like  the  tubes  with  rectangular  fins,  most  of  these 
tubes  displayed  enhancements  above  the  area  ratio,  as  shown 
by  the  two  Eo/ft r  columns  in  Table  II.  These  ratios  range 
from  0.94  to  2.95  under  vacuum  conditions,  and  from  1.52  to 
4.13  at  atmospheric  pressure.  Like  tube  number  15,  the 
small  values  of  0.94  and  1.52  for  the  "Wolverine”  tube 
resulted  from  the  large  amount  of  retained  condensate,  and 
by  the  very  thin  fins  (about  0.3  mm)  as  discussed  earlier. 
As  was  the  case  for  rectangular  fins,  the  highest  Eo/Ar  of 
2.95  and  4.13  were  shown  for  the  tube  (tube  number  36)  with 
a  relatively  small  area  ratio  of  1.37,  indicating  once  again 
the  importance  of  surface-tension  effects. 

3 .  ^e  Performance  of  " Parabolic”  Fins 

In  order  to  study  the  parabolic  fin  more  closely,  a 
tube  with  parabolic  fins  was  manufactured  (tube  number  38} 
with  the  same  fin  height,  spacing,  and  thickness  as  tube 
number  17  (tested  by  Georgiadis  [7]  and  this  investigator). 
The  performance  of  these  two  tubes  is  compared  in  Figure 
5.12  under  vacuum  conditions,  and  in  Figure  5.  13  at  atmos¬ 
pheric  pressure.  Clearly,  the  parabolic  fins  provide  a 
greater  enhancement  than  the  rectangular  fins.  The  reason 
for  this,  as  mentioned  earlier,  is  probably  the  continuous 
change  of  radius  of  curvature  (increasing  from  the  trough 
area  to  fin  tip)  for  the  parabolic  shape.  The  resulting 
pressure  gradient  along  the  fin  height  is  more  favorable 
(than  on  rectangular  fins),  which,  in  turn,  leads  to 
improved  drainage  from  the  fin  surface.  This  effect  was 
seen  most  clearly  on  tube  number  37.  Another  observation 
(that  contributes  to  improved  heat  transfer)  made  on  the 
tube  with  "parabolic”  fins  is  that,  although  completely 
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flooded,  the  fin  tips  (about  half  the  fin  height)  were 
visible  along  the  majority  of  the  fin  circumference  (except 
for  about  the  bottom  30  degrees) . 

As  discussed  in  section  E.1,  tube  number  37 
performed  poorer  than  the  spirally  finned  tube  (tube  number 
36)  under  vacuum  conditions.  The  reasons  for  this  are  not 
fully  understood;  so,  further  tests  of  these  fin  geometries 
(i.e.,  spiral  and  parabolic)  should  be  performed. 

P.  EPFPCT  OF  FIH-HETAI  THERMAL  COHDOCTIVITY  OH  STEAM-SIDE 

COEFFICIEHT 

To  investigate  the  effect  of  fin-metal  thermal  conduc¬ 
tivity  on  the  heat- transfer  performance,  a  stainless  steel 
tube  (tube  number  40}  was  manufactured  with  the  same  fin 
dimensions  as  the  "optimum"  copper  tube  (tube  number  6)  with 
rectangular  fins,  and  a  fin  height  of  1.0  mm.  Due  to  the 
low  thermal  conductivity  of  stainless  steel,  it  was  neces¬ 
sary  to  reduce  the  tube  wall  thickness  to  0.5  mm  to  prevent 
an  excessive  wall  thermal  resistance.  A  copper  tube  (tube 
number  39)  with  the  same  dimensions  as  the  stainless  steel 
tube  was  manufactured  for  comparison.  Also,  two  smooth 
tubes  having  an  outside  diameter  equal  to  the  root  diameter 
of  the  finned  tubes  (tubes  number  41  and  42}  were  manufac¬ 
tured;  the  first  tube  was  made  of  copper,  and  the  other  was 
made  of  stainless  steel.  The  results  for  data  runs  taken 
under  vacuum  conditions,  and  at  atmospheric  pressure  are 
shown  by  Figures  5.14  and  5.15  respectively. 

As  expected,  the  copper  finned  tube  provided  a  signifi¬ 
cant  enhancement,  while  the  stainless  steel  finned  tube 
actually  reduced  the  heat-transfer  performance  slightly  in 
comparison  to  the  smooth  tube  (tube  number  42)  .  The  thermal 
conductivity  of  stainless  steel  is  much  lower  than  that  of 
copper  (i.e.,  15  V/m.K  compared  to  385  W/m.K),  so  the 


Figure  5.14  comparison  of  Heat-Transfer  Performance  Between 
Stainless  Steel  and  Copper  Tubes  (Vac.  Runs) . 


Figure  5.15  Cosparison  of  Heat-Transfer  Perforaance  Between 
Stainless  Steel  and  Copper  Tubes  (&ta.  Runs) . 


thermal  resistance  of  the  stainless  steel  fins  is  signifi¬ 
cantly  higher,  resulting  in  a  poorer  fin  efficiency.  This 
is  combined  with  the  deleterious  effects  of  retained  conden¬ 
sate  (introduced  by  the  use  of  fins)  so  that  the  above 
factors  dominate  over  the  beneficial  thinning  of  the  conden¬ 
sate  film  between  fins.  The  net  effect,  then,  is  a  slight 
degradation  of  the  heat- transfer  performance  compared  to  a 
smooth  tube. 


G.  EFFECT  OF  STEAH  VEIOCITF 

As  discussed  in  Chapter  II,  most  models  of  film  ccnden- 
sation  on  finned  tubes  have  neglected  the  effect  of  vapor 
shear.  For  this  reason,  a  decision  was  made  to  perform  runs 
for  the  ’’optimum”  tube  (tube  number  6)  with  rectangular 
fins,  and  fin  height  of  1.0  mm,  at  nominal  steam  velocities 
of  2,  4,  6,  and  8  m/s.  Steam  velocity  through  the  test 
section  was  controlled  by  the  boiler  input  power,  and  since 
a  steam  velocity  of  1  m/s  was  the  maximum  obtainable  at 
atmospheric  pressure,  the  data  runs  were  performed  under 
vacuum  conditions  only.  The  results  of  data  taken  for  the 
above-mentioned  steam  velocities  are  shown  in  Figure  5.16, 

The  trend  shown  in  Figure  5.16  indicates  that  increased 
steam  velocity  (i.e.,  increased  vapor  shear)  leads  to 
enhancements  of  up  to  10  percent  as  steam  velocity  increases 
from  2  to  8  a/s.  This  is  to  be  expected  since  the  shear 
stress  at  the  liquid-vapor  interface  tends  to  pull  the 
condensate  toward  the  bottom  of  the  tube,  thereby  thinning 
the  condensate  film  at  the  top  part  of  the  tube. 

The  above-mentioned  data  were  compared  to  data  taken  by 
Tau  et  al.  [1]  for  steam  velocities  of  0.5,  0.7,  and 
1.1  m/s  at  atmospheric  pressure.  They  reported  enhancements 
of  up  to  40  percent  for  the  steam  velocity  of  1.1  m/s 
compared  to  the  0.5  m/s  velocity.  The  rather  large 
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difference  in  enhancements  reported  by  Yau  et  al.  compared 
to  those  of  this  investigation  was  not  expected.  It  should 
be  kept  in  mind,  however,  that  the  data  of  Yau  et  al.  were 
collected  at  atmospheric  pressure  (compared  to  vacuum  condi¬ 
tions  for  this  investigation),  where  the  increased  satura¬ 
tion  pressure  leads  to  a  lower  condensate  viscosity  (about 
45  percent  as  discussed  in  section  E. 1) .  This,  in  turn, 
allows  condensate  to  flow  more  freely  so  that  the  effect  of 
shear  stresses  acting  on  the  liquid-vapor  interface  would  be 
more  pronounced.  Another  factor  to  be  considered  is  the 
slightly  larger  fin  spacing  used  by  Yau  et  al.  (2.0  mm 
compared  to  1.5  mm  for  this  investigation).  If  the  fin 
spacing  is  very  small,  the  vapor  shear  may  not  be  "felt”  on 
the  sides  of  the  fins  or  on  the  fin  bases.  This  is  due  to 
the  fact  that  the  steam  flows  around  the  tube,  thus  leaving 
nearly  stagnant  steam  between  fins.  As  the  fin  spacing 
increases  though,  the  fin  sides  and  bases  will  be  more 
"accessible"  to  steam,  therefore,  giving  rise  to  the  benefi¬ 
cial  effects  of  vapor  shear.  Thus,  the  larger  enhancements 
of  Yau  et  al.  [1]  compared  to  this  investigation  may  be 
explained  by  the  two  phenomenon  just  described. 

An  attempt  was  made  to  correlate  the  data  using  a 
Fujii-type  [29]  equation  given  by: 

F.e,p0-^  (5.3) 
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Vg  =  steam  velocity 

=  steam-side  heat- tra nsf er  coefficient  based 
on  actual  surface  area  (including  fins) 

Rej.p  =  two-phase  Reynolds  number 

c,  d  =  constants  to  te  determined  experimentally 


The  constants  c  and  d  used  in  the  Fujii  equation  [29],  which 
was  developed  to  account  for  vapor  shear  on  smooth  tubes, 
were  0.96  and  0.2  respectively.  Values  of  c  =  1.302  and 
d  =  0.236,  however,  were  found  to  more  closely  fit  the  data 
of  this  investigation  after  a  least-squares  fit  was 
performed.  A  comparison  of  the  experimental  data  to  this 
curve  is  shown  in  Figure  5.17.  Although  a  fairly  close 
agreement  has  been  obtained,  the  orientation  of  data  clus¬ 
ters  does  not  match  the  slope  of  the  least-squares  line 
generated  by  globally  fitting  all  the  data.  It  is  possible 
that  the  errors  in  the  experimental  constants  used  in  the 
Sieder-Tate-type  equation  are  partly  responsible  for  this 
questionable  trend.  Georgiadis  [7]  showed  that  using  larger 
values  (from  the  "modified  Wilson  method")  of  the 
5)ieder-Tate  coefficient  (C)  tends  to  rotate  the  data  clus¬ 
ters  to  more  closely  match  the  experimental  correlation. 
Further,  a  Fujii-type  correlation  may  be  inappropriate  for 
finned  tubes  due  mainly  to  condensate  retention  that  occurs. 


note  that  the  slope  of  the  line  predicted  by  equation 
(5.3)  for  the  experimental  data  does  not  match  the  slope  of 
the  Fuji!  equation  mentioned  earlier,  but  instead,  nearly 
matches  the  slope  of  the  Husselt  equation  for  smooth  tubes 
gi7en  by: 
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which  can  be  rewritten  as: 
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(5.8) 


In  this  manner,  the  effect  of  vapor  shear,  at  least  on  the 
tube  tested  (tube  number  6)  ,  is  very  small  compared  to 
smooth  tubes  [7].  In  fact,  using  the  experimentally-found 
d-value  in  equation  (5.3),  the  steam-side  coefficient  and 
steam  velocity  are  related  as  follows: 


.0. 028 


(5.9) 


Once  again,  this  observation  shows  that  the  effect  of  steam 
velocity  on  the  steam-side  hea t-transf er  coefficient  is  very 
small . 

A  vertical  shift  of  the  least-squares-f it  compared  to 
the  Nusselt  theory  results  from  the  insreased  c-value  (1.30  2 
compared  to  0.728)  of  equation  (5.8).  In  fact,  the  ratio  of 
the  experimental  c-valie  to  the  coefficient  in  the  Nusselt 
equation  (i.e.  ,  1.302/0.  728  =  1.788)  represents  the  heat- 
transfer  enhancement  ii,  excess  of  the  area  increase  due  to 


finning.  Therefore,  the  s-value  depends  fin  geometry,  and 
surface  tension  of  the  condensate  etc.. 

In  view  of  the  above  discussion,  further  attempts  to 
modify  the  Nusselt  equation  to  include  vapor  shear,  surface- 
tension  effects,  and  fin  geometry  ate-  may  lead  to  an  accu¬ 
rate  correlation  for  the  steam-side  heat-transfer 
coefficient  for  finned  tubes. 

H.  EFFECT  OF  ZHTERHll  AND  EXTERNAL  ENHANCEMENTS  ON  THE 

OVERALL  COEFFICIENT 

Dp  to  this  point,  all  discussions  have  delt  with 
enhancements  to  the  steam-side  heat-transfer  coefficient. 
Even  with  large  steam-side  enhancements,  however,  a  large 
thermal  resistance  through  the  tube  wall,  or  on  the  water 
side  may  dominate,  and  completely  negate  the  steam-side 
enhancement.  For  this  reason,  a  decision  was  made  to  test 
the  "optimum”  tube  with  rectangular  fins  (fin  height  of 

I. 0  mm  -  tube  number  6),  and  the  corresponding  smooth  tube 
(tube  number  1)  ,  with  and  without  the  insert.  The  results 
of  these  data  runs,  taken  under  vacuum  conditions,  are  shown 
in  Figure  5.18.  At  a  water  velocity  of  about  3  m/s,  use  of 
the  insert  in  the  smooth  tube  resulted  in  a  40  percent 
increase  of  the  overall  coefficient,  as  did  the  use  of 
external  fins  without  the  insert.  By  combining  the  use  of 
fins  and  the  insert,  an  enhancement  of  nearly  155  percent 
was  attained.  It  must  be  noted,  however,  that  the  insert 
used  during  these  experiments  is  not  suited  for  use  in 
condensers  due  to  the  very  high  pressure  drop  it  creates. 
Instead,  Hebb  et  al-  [30]  showed  that  the  use  of  internal 
flutes,  fins,  or  circumferential  ribs  may  be  considered  to 
increase  the  water-side  coefficient  by  a  factor  of  up  to  2.0 
(compared  to  2.  5  obtained  by  the  insert  used  during  this 
investigation).  As  a  result,  it  appears  that  the  use  of 
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external  fins  togather  with  internal  enhancements,  could 
increase  the  overall  heat- transfer  coefficient  by  approxi¬ 
mately  100  percent. 

I.  DEYEIOPHERT  OF  A1  EZPEBIHEHTAL  COBBELATIOI 

Using  the  1985  model  of  Webb  et  al.  [25]  (discussed  in 
Chapter  II)  as  a  starting  point,  attempts  were  made  to 
correlate  the  experimental  data  obtained  for  tubes  with 
rectangularly-shaped  fins.  Time  limitations  prevented  an 
in-depth  effort,  but  a  few  simple  modifications  were  made. 
Before  discussing  these  modifications,  the  general  scheme  of 
the  Webb  et  al.  model  will  be  discussed.  A  more  detailed 
step-by-step  discussion  is  provided  in  Appendix  A. 

First,  a  value  for  the  parameter  ^  must  be  determined 
for  the  fin  geometry  in  question.  With  this  value  known, 
the  heat-transfer  coefficient  for  the  unfloodei  fin  surfaces 
can  be  obtained  using  equation  (2.30).  Values  for  the  heat- 
transfer  coefficients  of  the  unflooded  tube  area  between 
fins,  and  for  the  flooded  region  are  determined  next  using 
equations  (2.32)  and  (2.33)  respectively.  The  above- 
mentioned  heat- transfer  coefficients  are  then  substituted 
into  equation  (2.31)  to  obtain  the  average  heat-transfer 
coefficient  for  the  finned  tube  in  question. 

A  careful  review  of  the  Webb  et  al.  model  revealed  two 
deficiencies  worth  noting.  First,  in  their  two-dimensional 
analysis  of  the  flooded  tube  region,  the  heat-transfer  coef¬ 
ficient  for  the  flooded  root  area  between  fins  was  applied 
to  the  fin  tips.  Figure  2.4  shows  that  a  constant,  rela¬ 
tively  thick  condensate  film  is  assumed  which,  when  applied 
to  the  fin  tips,  results  in  an  underprediction  of  the  heat 
transfer.  Use  of  the  Nusselt  equation  for  smooth  tubes 
would  be  more  appropriate  for  fin  tips  in  this  region. 
Additionally,  the  effects  of  convection  in  this  region  are 


ignored;  however,  to  include  these  effects  would  make  a 
solution  extremely  complex. 

The  second  deficiency,  probably  valid  only  for 
rectangularly-shaped  fins,  involves  the  assumed  condensate 
film  shape  between  fins  on  the  unflooded  tube  surface. 
Figure  2.4  shows  a  depression  in  the  film  next  to  the  fin 
base,  with  an  increasing  condensate  thickness  away  from  it. 
Condensate  film  shapes  observed  during  this  investigation 
continued  to  decrease  away  from  the  fin  base  as  shown  by 
Figure  2.1,  probably  due  to  the  thinning  effect  of  surface- 
tension  forces.  Based  on  the  assumption  of  a  relatively 
thick  condensate  film,  predictions  using  the  nebb  et  al. 
model  will  underpredict  heat  transfer  from  this  part  of  the 
tube. 

To  use  the  Webb  et  al.  model  for  rectangular  fins,  a 
value  of  -0.9  was  chosen  for  5  .  A  value  of  -1.0,  which 
corresponds  to  a  straight  line  (i.e.,  straight  side  of  the 
rectangular  fin)  would  seem  more  appropriate;  however,  this 
value  would  lead  to  zero  heat  transfer  from  the  fin  sides 
using  equation  (A.1).  Webb  et  al.  used  values  around  -0.7 
to  -0.8  for  their  slightly  trapezoidal  fins;  so  the  value  of 
-0.9  was  assumed  to  be  reasonable.  Using  this  value  along 
with  a  4)  value  of  5.0  (used  by  Webb  et  al.)  resulted  in 
significant  underprediction  of  the  experimental  data 
obtained  by  this  investigator  and  by  Georgiadis  [7],  as 
shown  in  Figure  5.19. 

In  an  attempt  to  improve  the  predictions,  values  for  the 
41  ratios  were  calculated  using  the  paralled  heat  flow  paths 
prnposed  by  Owen  et  al.  [13]  discussed  earlier  in 
Chapter  II.  These  values  are  dependent  on  fin  spacing, 
and  should  decrease  with  increasing  fin  spacing.  This  trend 
was  observed,  and  for  the  set  of  tubes  with  both  fir.  height 
and  fin  thickness  of  1.0  mm  (s=  0.5,  1.0,  1.5,  2.0,  and 
4.0  mm  -  tubes  number  4  to  8),  resulted  in  4>  values  that 
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ranged  from  1.7  (s  »  0.5  nm)  to  1.2  (s  =  4.0  mm).  Using  the 
Febb  et  al.  model  once  again,  but  with  a  constant  average 
value  of  1.5,  gave  the  results  shown  in  Figure  5.29  for  the 
above-mentioned  tubes  (tubes  number  4  to  8) .  This  Figure 
also  shows  a  plot  of  the  Hebb  et  al.  model  using  the  vari¬ 
able  (i.e.,  1.7  to  1.2)  which  made  almost  no  difference. 

Cnee  again,  the  experimental  data,  also  shown  in  Figure  5.20 
were  significantly  underpredicted. 

In  order  to  correct  for  the  assumed  thick  condensate 
layer  between  fins,  the  work  of  Fujii  et  al.  [31]  was  used. 
Their  paper  delt  with  smooth  tube  enhancements  from  coiling 
a  wire  around  the  tube  circumference,  and  provided  an  equa¬ 
tion  to  account  for  the  thinning  effect  between  wire  wraps 
as  shown  below: 
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(5.  10) 


where 


=  wire  diameter  (taken  as  the  fin  height 
for  use  in  the  Febb  et  al.  model) 


(5.11) 


(5.12) 


Adding  this  correction  factor  to  the  heat-transfer  coeffi¬ 
cient  for  the  unflooded  root  area  between  fins  made  a 
significant  Improvement  for  the  group  of  tubes  mentioned 


earlier  (tubes  number  4  to  8)  as  shown  by  Figure  5.20. 
Nonetheless,  experimental  data  were  still  underpredicted. 

In  summary,  it  should  be  kept  in  mind  that  this  model 
was  not  developed  for  rectangular  fins.  The  most  signifi¬ 
cant  discrepancy,  however,  appears  to  be  on  the  prediction 
of  heat  transfer  through  the  flooded  area.  This  model 
predicts  only  minimal  heat  transfer  through  the  flooded 
portion  of  the  tube,  whereas  Panniarachchi  et  al.  [2] 
reported  enhancements  at  least  equal  to  the  area  ratio  for 
completely  flooded  tubes.  Additionally,  there  may  be  other 
phenomena  which  have  not  been  taken  into  account  (such  as 
intense  convection  in  the  condensate)  which,  if  included, 
may  improve  the  predictions.  Nonetheless,  this  model  shows 
promise  and  should  be  pursued  further. 


VI.  COHClPSlOHS  MS  SBCOaHEHDATIOHS 


coiciasiovs 

.  The  use  of  fins  may  lead  to  significant  enhancements 
of  the  steam-side  heat -transfer  coefficient.  Under 
vacuum  conditions,  enhancements  as  high  as  4.0  (tube 
number  36)  were  realized,  vhile  enhancements  of  up  to 
6.9  (tube  number  37)  were  realized  at  atmospheric 
pressure.  These  enhancements,  which  were  up  to  three 
times  greater  than  the  area  ratio  (finned  tube 
area  /  smooth  tube  area)  are  due  to  the  thinning 

effect  of  surface-tension  forces. 

.  For  rectangularly-shaped  fins  with  fin  heights  of  0.5 
and  1.5  mm,  the  2.0  mm  fin  spacing  appears  to  be  the 
optimum  among  the  fin  spacings  tested  (1.0,  1.5,  2.0, 
and  4.0  mm).  This  was  the  case  for  both  vacuum  and 
atmospheric  conditions.  For  similar  tubes  with  fin 
heights  of  1.0  and  2.0  mm,  Georgiadis  *7]  found  the 
optimum  spacing  to  be  1.5  mm. 

.  The  use  of  fins  on  a  stainless  steel  (low- thermal- 
conductivity-  metal)  tube  degraded  the  heat-transfer 
performance  in  comparison  to  its  corresponding  smooth 
tube.  In  this  case,  the  poor  fin  effeciency  combined 
with  the  deleterious  effects  of  retained  condensate 
dominated  over  the  beneficial  thinning  of  the  conden¬ 
sate  film  between  fins,  resulting  in  a  net  reduction 
of  the  heat- transfer  performance. 

.  Fin  geometry  affects  the  extent  to  which  the  surface- 
tension  forces  thin  the  condensate  film  on  the  fin 
surfaces.  Parabolically-shaped  fins  were  shown  to 


oaximize  the  thinning  effect  of  surface-tension 
forces. 

.  Increasing  steam  velocity  (i.e.,  vapor  shear)  from  2 
to  8  m/s  increased  the  steam-side  coefficient  by 
about  10  percent  (under  vacuum  conditions) ,  in 
contrast  to  7au  et  al.  [ 1 ]»  who  showed  a  40  percent 
increase  when  steam  velocity  was  increased  from  0.5 
to  1.1  m/s  (at  atmospheric  pressure). 

.  For  the  "optimum”  tube  with  rectangular  fins  and  fin 
height  of  1.0  mm  (tube  number  6)  ,  the  overall  heat- 
transfer  coefficient  without  the  insert  showed  an 
enhancement  of  40  percent  over  the  corresponding 
smooth  tube.  With  the  insert  in  place,  the  finned 

tube  showed  overall  enhancements  of  about  155 

percent. 

.  Host  theoretical  models  developed  thus  far  to  predict 
the  heat-transfer  performance  of  steam  condensation 
on  horizontal  finned  tubes  neglect  one  or  more 
phenomena  of  major  importance.  Thus,  there  are  no 
models  currently  available  to  accurately  predict  the 
performance  of  finned  tubes  used  to  condense  steam. 

BECOHBENDATIOHS 

.  Hake  further  refinements  to  the  Webb  et  al.  [25] 
model  in  an  effort  to  predict  the  experimental  data 
within  a  reasonable  accuracy. 

.  Manufacture  and  test  additional  finned  tubes  with 
parabolic  and  spiral  fin  shapes  with  different  dimen¬ 
sions  (i. e.,  fin  spacing,  height,  and  thickness).  In 
this  manner,  determine  the  optimum  dimensions  for 
each  fin  shape. 


Manufacture  and  test  a  series  of  tubes  with  varying 
tube  diameter  to  study  the  effect  of  diameter  on  the 
heat-transfer  performance.  For  this  purpose,  select 
the  optimum,  rectangular-shaped  fins  (s  =  1.5  mm, 
e  =  1.0  mm,  and  t  =  1.0  mm). 

Collect  data  on  finned  tubes  fitted  with  drainage 
strips  made  out  of  both  solid  and  porous  metal  to 
determine  the  enhancements  they  provide. 

Continue  tests  on  finned  tubes  with  tube-metal 
thermal  conductivities  between  those  of  copper  and 
stainless  steel,  such  as  copper-nickel  alloys  and 
aluminum. 

Conduct  data  runs  at  atmospheric  pressure  with  steam 
velocities  of  0.5,  0.7,  and  1.1  m/s,  using  tube 
number  6  to  compare  with  the  data  of  Tau  et  al.  [1]. 

Perform  data  runs  using  different  fluids  with  widely 
varying  fluid  properties  (i.e.,  surface  tension, 
viscosity,  etc.)  to  further  study  these  effects,  and 
to  compare  with  existing  theoretical  models  and 
experimental  correlations. 


APPENDIX  A 

PROCEDOBE  FOR  OSE  OP  THE  BEBB  ET  AL.  [25]  HODEl 


As  proposed  by  Webb  et  al.,  solving  for  the  steam-side 
heat-trar.sfer  coefficient  proceedes  as  follows: 

1.  First,  a  value  fore  must  be  determined.  This 
parameter,  as  shown  Ly  figure  2.4,  is  a  measure  of 
the  fin  aspect  ratio  (e/t) ,  and  it  must  be  deter¬ 
mined  using  an  iterative  method.  It  is  a  function 
of  fin  geometry  only,  and  should  match  the  fin  side 
profile  as  closely  as  possible.  Webb  et  al.  [25] 
provide  details  of  this  method,  and  no  further 
discussion  is  provided  on  this  since  a  true  ^-value 
does  not  exist,  as  defined  by  Adamek  [27],  for  a 
rectangular-shaped  fin. 

2.  With  the  value  of  known,  determine  the  heat- 
transfer  coefficient  for  the  fin  sides  in  the 
unflooded  portion  of  the  tube.  This  value  is 
obtained  using  the  following  equation  [25,  27]: 


hj,  =  2.U4 
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3.  Determine  the  heat-transfer  coefficient  for  the 
unflooded  root  area  between  fins  using  the  Nusselt 
equation  written  in  terms  of  the  film  Reynolds 
number.  For  this  purpose,  use  the  following  itera¬ 
tive  method: 


a)  Calculate  the  condensation  rate  (m)  using: 


r 


(A.2J 


Assume: 


=  m 


(A. 3) 


c)  calculate  the  film  Reynolds  number  by: 


Re  =  _ £ _ 

^  Pf  fCs  +  t^')  -  t>,l 


(A. 4) 


d)  Calculate  h  using: 
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h  =  1.514  JiL- — 
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(A. 5) 


e)  Calculate  the  base  heat  flow  by; 


^  =  »V  AT 


(A. 6) 


f)  Solve  for  m  using: 


•  '  'r* 
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(A. 7) 


g)  Repeat  steps  c  thru  f  until  convergence  is 
obtained. 
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Detecmination  of  the  heat- transfer  coefficient  in 
the  flooded  part  of  the  tube  requires  that  a 
computer  solution  be  used  to  solve  the  two- 
dimensional  conduction  problem  of  heat  flow  through 
the  fins  and  retained  condensate.  This  solution 
provides  a  ratio  ((|>)  of  actual  heat  flow  into  the 
coolant  divided  by  the  heat  flow  that  would  result 
from  fins  of  zero  thickness.  The  heat-transfer 
coefficient  is  then  determined  by: 

\  -  '^-L  (A. 8) 

f* 


The  fin  efficiency  is  determined  next  using  the 
following  equation: 

^  _  tanh  fm  (e  +  t/2>] 

m  (a  +  t/?.>  (A*9) 

where 

_  -1/2 
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(A- 10) 

.  ™  . 

It  should  he  kept  in  mind  that  this  fin  efficiency 
is  only  valid  for  "thin”  fins.  Nonetheless,  its  use 
is  justified  since  there  are  no  equations  available 
specifically  for  "thick"  fins. 
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APPENDIX  B 

HODIPIED  NIISOM  METHOD 


The  •’modified  Nilson  method"  calculates  the  leading 
coefficient  foe  the  Sieder-Tate  equation.  This  method 
assumes  a  form  of  correlation  for  both  the  water-side  and 
steam-side  heat-transfer  coefficients,  each  of  which 
contains  a  coefficient  to  be  determined  by  iteration.  In 
the  past,  a  Musselt-type  equation  was  used  for  the  steam- 
side  heat-transfer  coefficient  as  shown  below  [7]: 
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Equation  (B.1)  results  in  a  6  value  of  0.655  for  conditions 
where  no  vapor  shear  is  present.  The  presence  of  vapor 
shear  generally  results  in  a  higher  value,  which  must  be 
determined  iteratively.  The  use  of  this  equation  should, 
therefore,  be  restricted  to  relatively  low  vapor  shear 
conditions  (i.  e.  ,  vapor  velocity  less  than  1.0  m/s)  .  Note 
that  a  Nusselt-type  equation  does  not  represent  the  "correct 
form,"  as  it  contains  no  terms  involving  vapor  shear.  To 
avoid  this  deficiency,  a  correlation  developed  by  Fujii  and 
Honda  [29]  which  takes  the  effect  of  vapor  shear  into 
account  was  used,  and  is  shown  below: 


1/5 


Rewritiny  equation  (B.2)  to  express  h  as  a  function  of  heat 
flux  and  vapor  velocity  results  in: 


ho  »  B 
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(B.3) 


The  constant  g  in  equation  (B.3)  is  to  be  determined  by 
iteration.  A  Sieuer-Tate-type  equation  is  used  for  the 
water-side  hea t- transfer  coefficient,  as  shown  below: 
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Equations  (B.3)  and  (B.4)  are  then  substituted  into  the 
equation  for  the  overall  heat-transfer  resistance  given 
below : 
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A  linear  equation  used  to  generate  the  Wilson  plot  is 
obtained,  and  shown  below: 
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The  parameters  X  and  T  are  determined  from  fluid-property 
values  and  heat-flux  measurements  taken  during  data  runs 
using  an  uninstrumented  smooth  tube.  Iterations  between  the 
coefficient  C,  and  the  coefficient  B  ,  are  continued  until 
convergence  of  the  coefficients  (between  two  successive 
iterations)  within  0.1  percent  occurs.  The  slope  of  the 
Wilson  plot  generated  is  the  reciprocal  of  the  desired 
Sieder-Tate  coefficient. 
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APPENDIX  C 
LISTIN6  OF  EAV  DATA 


The  following  pages  contain  raw  data  obtained  for  tubes 
nuober  26  thru  44  under  vacuum  conditions  and  at  atmospheric 
pressure.  Raw  data  for  tubes  number  1  thru  25  are  presented 
by  Georgiadis  [  7  ]. 
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iPPENMI  5 

OHCBST&INTT  &MALTSIS 


There  is  always  an  uncertainty  associated  with  any  neas- 
urement  which  is  dependent  on  the  measuring-device  accuracy, 
as  well  as  on  the  operator’s  experience.  Numerical  data 
collected  during  this  thesis  effort  were  used  together  with 
theoretical  formulations,  so  final  values  of  the  steam-side 
heat-transfer  coefficient  may  be  distorted  due  to  uncer¬ 
tainty  propagation  during  calculations.  In  cases  where  the 
final  results  show  large  uncertainties,  it  may  be  unwise  to 
draw  any  conclusions.  Instead,  the  apparatus  and/or  the 
measuring  techniques  may  have  to  be  modified. 

Estimates  of  uncertainties  for  this  investigation  were 
obtained  using  program  ''UNA6”  which  is  listed  at  the  end  of 
this  appendix,  along  with  selected  uncertainty  evaluations. 
This  program  determines  the  uncertainties  using  an  equation 
proposed  by  Kline  and  McClintok  [32]  shown  below: 


where 

is  the  uncertainty  of  the  desired  dependent 
variable 

,  X2 ,  ....,  x^  are  the  measured  (independent) 
variables 

Wj^  ,  w^,  ....,  w^  are  the  uncertainties  in  the  measured 
variables 


'lAr 


!  I'.i'l*  -"’LI 

MI’O!  RlVIZZ:  :  -^-cr 

ic::o' 

•o^o  COM  /Cc'  ■:</') 

UUO  DIM  E(4> 

lOSu  DftCa  D.  lOOfHSOb;  ,:'Or27. 54363,-/ 67:40.3295.78025595. 8’ 

1060  DftTA  -9247486589. 6. 97b33£+ll .-2. 66'92E*n. 3. 94078E+14 
IC7'J  <£9D  <’'(•> 

1080  ,=R!MT  USING  "I  OX  .  ••"DATA  PQR  thE  UNCERTAINTY  ANALYSIS: . 

1090  PRIN"^ 

1100  PRINTER  I.S  701 
’110  BEEP 

M:■n  INPUT  "ENTER  PILE  NAIlE".Pile5 

1130  ARINf  USING  •■10X."”Fiie  Nane  :  .  1  2A"  :Fi  leS 

I  1  /lO  BEE-* 

1150  INPUT  "ENTER  DATA  SET  NUMBER  FOR  UNCERTAINTY  ANAL YSTS" . Ids 
1160  BEEP 

1-7Q  INPUT  "ENTER  PRESSURE  C3MDITIGN  »0-V.1-A>  ".Prc 
1171  Prc-Pfc*'' 

MOO  ASSIGN  iPi'e  TO  PMab 

1190  ENTER  AF  ;  l6  ;  I f g  .  I'm 

I19'  IF  It'3*U  then  Enter  •File:Dd 

1192  IF  Ifg*l  THEN  ENTER  »F i le ; Do . Dd .Dd  ) 

1200  FOR  1*1  TO  la; 

■  2 ' 0  ENTER  .“F  1  le  :  Bvo  I .  Banp  .  Utran  .£  tp  .£<  ■•  >  .^n.  Tci .  "co  .PSg .  P water 
1220  NEXT  I 
1230  Enf-E(O) 

’240  IP  Prc-1  THEN 
1 250  BEEP 

'250  TJcjTfi]-  ijSifiG  '■  1  OX  ,  ""Pressure  Condition:  'vacuum  <11  <Pa) . 

i:7G  E'.SE 

1250  PRINT  'j3I,''IG  "IOY  .""'-ressure  'lonaition.  Atnosoneric  ('01  kPa)'"’ 

•29','  £N0  IF 

ISOO  PAINTER  IS  1 
I? 10  BEEP 

I  320  PRINT  USI.NG  "4X. ""Select  firie  -ua  i  i  tyoe:  0»rriicK  I'tnm . 

•??D  INPU*  Itt 
'.‘vij  IE  itf'i  Then 
i3Si';  Do*.  01 905 
1 360  D ! • . U 1 27 
1370  ELSE 
II'BO  Do-.  01 45 
I39ij  D1-.OI35 
1 41.10  END  IF 
lain  PRINTER  19  1 
1 420  PEEP 

i43u  PRINT  L'SItiG  "4X  ,  ""ijelec’  natenal  code: . 

1640  PRINT  USING  "4X,""0-Cu,  '••SS . 

1450  INPUT  Inc 
1460  IF  Im.c-O  THEM 
1470  k'c-335 

’471  Orc-10 
1400  Cl-. 07’ 

14AL'  Ac-0 
1500  ELSE 
'510  Xc-15 

1511  Dlrc-1 
1520  Cl -.0686 
1530  Ac-0 
'540  END  IF 
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1^50  KRIMTE9  [3  701 
15(30  Ts*^  1tTvsv<  Enf  1 

[570  PRI.-lf  USIMG  "  I  OX  . ""Gteaw  7»moi?ratufa  •  .  i'D . 3D . <Deg  C> : 

I5j0  ORIfir  'jSINO  " ’  OX  . "-.Ucer  Fiou  Rate  O:)  -  . 30 . ID"  :Fr» 

1 5‘ifO  Dtc  i  •  .  0113 
1500  Dtco«.005 

1510  BEEP 
1 5Z0  Demt  *1 . 0E~6 

1630  0ts*6(JR<  <  <C<  I  >*2‘>'(2^-Emr  ?>•£."> f  4>-tmt  '  3> ’Denf  >  '2 > 

1540  T-<Tci+Tco)/2  I  FILM  'EMPESPnjRE 
:550»  oncer IfilNTY  IN  'HE  CQGlING  HrtTER 
1550  Drho-.5  !  ERROR  IN  UATEP  DEM'''nY 
1570  Dmt-.008  !  ERROR  IN  MRS'3  F..Q'h  RHIE 

1630  Rho-FNRho(T)  !  WATER  DEHSITY 

1530  Mf  - 1  .04805E-2*5.o0932E- -n  •  MASS  FLOW  RATE  QF  COOLING  WATER 
1  700*  CORRECT  MF  FOP  THE  TEMPERATURE  ErFECi" 

1710  Mf-MT  •<  1 .0365-!  .OSEAuE-a-Tci^S.I-SjE-o'Y.;^  2>  7  9354  la 
’720  Ai-(PI*0i- 2)/4  !  TUBE  INSIDE  CROSS  SECTION  AREA 
!7C0  Odi*. 000025 

!?40  'DaiJTPI^Odi  2)/4  1  ERROR  OF  INSIDE  TUBE  CROSS  AREA 
1750!  COMPUTE  THE  WATER  OELGCirr 

17F0  '/ij-Mf /(Rho*Ai )  '  WATER  OELOCITY  , 

1  770  PRINT  USING  "  I  OX  .-'Uatef  '/eiocitv  ■  "".2. 00.""  Tn/s)"'"'*;'/ 

KTSU!  CORRECT  OUTLET  ^^ATER  'E'^P .  ^CR  THE  ‘iIXING  CHAMBER  EFFECT 
1790  IF  Inn-1  OR  Inn-5  taen  'co>  T.ro- .  C'0a*«>y  3 
li-iOO  IF  r.-'n-O  then  ‘co- 'co-<  - .  00  1  is*  .  jO i -'I'tf  3) 

ISIO  T-<Tci*-Tco»-.5  '  FILM  TEmPF^ATijRE 

1330!  COMPUTE  THE  ERROR  IN  wATER  Of.JCITY 

1830  Owij-U.j-SOR<  <Dnt /Mf )  3*<  Drno/Pho*  2'KDai/Ai>  2) 

iSaO!  UNCERIAINTY  IN  THE  REYNOLDS  NUMBER 

1350  Mtj-fNMu<T)  !  WATER  VISCOSI'Y 

Omu-H.c-b  !  ERROR  0.'^  HATER  VISCOSITY 
1870  Re- ( Rnn-Vii)»D  i  > /Mw 

1  S30  Dre-Re-SQR ( < 0 r no/R5o )  '2 * •  O'/ij/ Vij )  2*!Ddi/Di)  2*iO")u7Mii()  2> 

1390!  UNCERTAINTY  IN  THE  HEAT  TRANSFERRED 
1900  Cp'u-FNCpij<  T ) 

!910  0-Nf  •  <  Tco- Tc i  >«Cpiu 
l'?3(i  Ocpia'O 

1330  0q-(3-‘>GR<  <  Ont /Mf  )  <  C  tco/ <  Tco- Tci )  •  >  2*<  <  Otci/(  Tco-’’c  1 1  •  >  2*<Dcpiii/'CpuJ  '2 

) 

1940!  UNCERTAINTY  IN  'HE  -EAT  FLUX 
!950  01 -.0005  !  ERROR  IN  TijsE  LENGHT 

1960  Ddo*. 000025 

'970  L-. 13335  !  CONDENSING  'URE  LENGTH 

1980  0o-0/<PI-Do-'.)  !  HEAT  PLUX 

1990  PRINT  USING  "  ' OX  . ""Heat  P'la*  -  .  Z  . 3DE . «W/ni  2> . : 

Qp 

3000  Dap- < Op/P  1 1 -SOR < < Dq/0 •  3*iDao/Do>  2**D1/L*  2> 

2010  Lmtd'i  Tco-TciJ/LQGi  <  T's-T.:!  //(  > 

3020  'Jo-Up/Liatd  !  OVESALi.  nEHT  TRANSFER  COEF. 

2030  A1-0ts-<Tci-Tco)/<<Ts-Tci)-<T-;-Tco>-L0G<<Ts-Tci)/(T3-Tco>(> 

3UaO  A2-0t>ti/<<T5-Tcii-'-UGT(T5-Tci)7(Ts-Tco>>) 

2050  Aj-Otco/iiTs-Tcol-cOGfiT-s-Tcil/iTs-TcoJi* 

2050  Olmrd-L.ntd-SOR* At  2*t*2  2+A3  2) 

2070  0uo-Uo-S0R<  <  Ddp/iJp  >  2*<  D  Imtd/Lmtd* '2> 

2050  M-Mu 

3090  T1-(T*37?.i5)/273.:5 

3100  Lw-FNKuiTl) 

3!;0  Ac-0.  !  INTERSCEPT  PROM  SIECcR  PROGRAM 


*«  DTIC  does  not 
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2120  LI  *.060325  1  LENGTH  Or  iJNEIfJNEr-  LEFT  PflPT  DF  THEE 
2130  L2-. 034325  !  lZ.NGTH  OF  iNFItmED  RIGHT  PORT  .3F  tiJEE 
2140  01*. 01905 

2!50  02-.0:5875 

2160  Pr*Cpu*Mii)/Ku 
2170  fluu*FNI1uu<  T ) 

21E10!  LiNCESTAINTY  OF  INSIDE  MEA7-TRANSFE.9  CQEFF. 

2190  Cf-1. 

2200  Hi-(Ku/Di)»<Ci-Ke'  .3*Pr  .333*Ct‘-Ac) 

2210  Dt i*O/<PI*0i “(L^L 1 “Fal *02*662  > “Hi ) 

2220  Cf c* < Nuu/FNMuui r*Dt : ) >  ' .  1 4 

2230  IF  PEtSMCfc-Cr>/Cfc>>  .01  THEN 

2240  Cf-<Cf *Cfc>*.5 

2250  GOTO  2200 

2160  END  IF 

2270  P1-PI»<Di*D1 > 

2230  Bl-<D1-Ln  )*PI*<Di*Oi >*.5 
2290  H1-(Hi*pi/(Kc*31))  .5 

2300  P2-PI«<0i+02) 

2310  B2-<D2-Di)*p:*<0i*D2>*.5 

2320  H2-<Hi*P2/i'<c*P2))  .5 

2330  Fe’-FNTanhifU-LI  i/^H’-LH 

2540  Fe2*^NTanh <H2*l2 > / < M2*L2 >  I 

2350  0tc*Q/<PI*Di*<L+Li*F6l  ♦■_2*'  62)*Hi  > 

2260  IF  ABS((Dtc-Oti)/Dtc)>.01  THEN  2200 

2370  Dir w. 0010  1  ERROR  IN  WATER  THERMAL  CONDUCTIVITY 

2330  Dci-.0a05  1  E.RROR  IN  SIEDER-tate  COEFFICIE.NT 

2390  Dor*. 05  !  ERROR  IN  PRANDTi.  NUMBER 

2400  Dcf*8.E-6 

2410  A4. , lA.Ocf /Cf 

2420  D*n*hi*SGR<  (Olfiii/Kw)  2*<0di/Di)  2*<  .3*Dr«/R6>  2*<  .332*0pr/Pr>  2*<Dcj/Cj  >  2* 
44  > 

2".50?  UNCERTAINTY  OF  OUTSIDE  HEAT- ■TRANSFER  COEFR. 

2440  Rw*Do*lOG(Do/Di )/<2*kc )  '  HALL  RESISTANCE 

2450  Ho*  ' /<  <  I /Uo ) - <  0o*L/ <  0 1  * < L**- 1  *6 6l  ♦L2*Fe2>*Hi  >  > -Ru) 

2460  DrwR«*S0R<  <Dao/Do )  ‘  2* <  Dk c/K  c » '■2*<D<Jo/'< Do*LOG< Do/Di  > ) ) '  2*(Ddi/<Dj*LOG(Do/D 
1 ' > )  2) 

2470  A5* 1 /Uo-Rw- < Do/ < D i “Hi J ) 

24S0  A6*Dao/<Uo- 2*A5) 

2490  A7*Cirij/H5 

2500  Afl*<  (Do/<Di*Hj  »> *<961/111 )  )/A5 
2510  PRINT 

2520  Dho-Ho*S0R<A6  2*A7'2*A3 '2) 

2530?  CALCULATE  THE  %  UNCERTAINTY  IN  Ho 
25^0  PrHo*06f>*  I  00/Ho 

2550!  CALCULATE  THE  7.  UNCERTAI.NTY  IN  REYNOLDS  NUMBER 
2560  Prre-Dre*1 00/Re 

25  70?  CALCULATE  THE  /.  UNCERTAINTY  IN  MASS  FLOW  RATE 
2530  PrnT "Dm’ *1 00/Mf 

2590?  calculate  THE  UNCERTAINTY  IN  hEAT  TRANSFER 
2600  Prqp*DaD*1 00/Gp 

2610?  CALCUL.ATE  'hE  UNCERTAINTY 'IN  LMTD 
2620  Pr  Imtd'D  intd*1  OO/'Lmtd 
2630?  CALCULATE  THE  /.  UNCERTAINTY  IN  Rw 
2640  Prru'Dru*  1 00/Riii 

2650?  CALCULATE  ThE  7.  UNCERTAINTY  IN  OVERALL  HEAT  TRANSFER  CCEF . 

2660  Pruc'*0iio*  1 00/Uo 

2670?  calculate  ThE  7.  UNCERTAINTY  IN  INSIDE  HEAT  TRANSFER  COEPF. 

2680  Pr6i*Dhi*100/H; 

2690  PRINT 

2700  PRINT  USING  "lOX. "•'UNCERTAINTY  ANALYSIS: . 
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2710 

PRINT 

2711 

PRINT 

USING  "lOX.""  UARIfiBLE 

percent 

UNCERTAINTY . 

2712 

PRINT 

2720 

PRINT 

USING  “I  OX.  •"•Mass  Flow  Rate.  Md 

•  ••• 

.2D.2D." 

iPrtnf 

2730 

PRINT 

USING  ••!  OX  .  ••••Reynoios  Ntunoer.  Re 

•e  •« 

.3D. 2D." 

tPrre 

2740 

PRINT 

USING  ••  I  OX.  ••••Heat  Flax.  A 

-.2. 2D. 

•;Prqp 

2750 

PRINT 

USING  ••10X. ••••Log-Mean-Ten  Diff.  LMTD 

-.2. 2D" 

:Pr Iwtd 

2  760 

PRINT 

USING  •'I  OX.  •••••Hall  Resistance.  Ru 

•.DD.2D. 

" :Prro 

2770 

PRINT 

USING  ••lOX.^’^'iTverall  H.T.C..  Uo 

*• 

•.DD.2D. 

" :Prao 

2700 

.PRINT 

USING  ••’OX. ■•••'Hater-Side  H.T.C..  Hi 

*••1 

.30. 2D." 

;Prhi 

2790 

PRINT 

USING  ••lOX.""Steam-Siae  H.T.C..  Ho 

•  •  •• 

.3D. 2D." 

:Prno 

2900 

END 

2810 

DEr  FNMuuO) 

2920 

A-247 

8/<T*l33. 15) 

2930 

Muo-Z 

.4E'5*10‘A 

2940  ;?ETURN  fluw 
2950  FNF.ND 
2iS0  DE‘  ^NrannOO 
2870  P-EXP<X> 

2880  (3-EXP<-X> 

2830  Tanr.-<P-0)/<P+0) 

2300  RETURN  Tann 
2910  EMEND 
2990  DEF  riVKuC?) 

2930  (<(j«-.32247  +  ti.<2.3395-ti*<  ;  .8007-TI*<  .52577-  . 07344-n  >  >  > 

2940  return  Ku 
2950  "MEND 
2960  DEF  "Nnw<T) 

2370  M-247.3/< T*I33. 15> 

2380  ‘t(j-2.4E-5*I0  A 

2990  return  tiu 
3000  FNE'ID 
3u10  OEF  "MRho<r) 

3020  860-999. 52946*T.<  ,  Jl2£9-''*<5.4825I3E-3-’’*!  .234147E-5> ) 

3030  RETURN  Rno 
3u40  FMEfiD 
3050  DEF  FNCou<T) 

3060  Cow- <4. 21 i20858-T-<2.26326E-?-T«<4.4236ft-5»2.?!429E-7-T) ) )-lOOO 

3070  RETURN  Cow 

3080  FMENO 

3050  DE'  rNTvSy<Eir.O 

31U0  COM  /Cc/  Cl  7) 

3110  T-C<C) 

3120  FOR  i-i  70  7 
3130  T-T+Cdi-Etnf  I 
3140  NEXT  ! 

3150  return  T 
3160  FNENO 
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Copy  available  to  DTIC  does  nc 
permit  fully  legible  reproduction 


FOR  the  iJNCERTAIMTY 

File  ^Jane : 

Pr-j-Saure  ('.oridition: 

Sttran  TeiTiDerature 
l^atar  Fiou  Rate  ('/> 

U'lter  'Velocity 
Heat  Flux 


UNtERlAiNTY  ANALYSIS: 
PARIAELE 

Mass  FIou)  Rate.  Md 
R'?yriolos  NiiPiber,  Re 
Heat  Flu-x.  q 
Lig-Mean-Teri  Ciff.  LMTD 
Hail  Resistance.  Ru 
Q''erai  1  H .  T  .  C .  ,  !Jo 
Water-Side  H.T.C..  Hi 
Steam-Side  H.T.C.,  no 


riQ'A  r.jR  7h£  uncertainty 

“ ! i e  Name ; 
f^'ressure  Condition: 

Steam  Temoeratiire 
Hater  Flow  Kate 
Hater  '•/‘elociry 
Heat  Flux 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Hass  Flow  Rate.  M(3 
ReyTiolds  Number.  Re 
Heat  Flux,  q 
Log-Mean-Tem  Oiff.  INTO 
Wall  Resistance.  Rw 
0  >erai 1  H . T .C. .  No 
Water-Side  H.T.C..  Hi 
Stean-Side  H.T.C.,  Ho 


ANALYSIS: 


F32U2G 

Uacuum 

<  ?  t  ^l^a) 

A.3.A7 

<  Deg  C ) 

20.00 

i  .  16 

i  m/s ) 

-  2.niE+05  <W/r,'2) 


PERCENT  UNCEPIAINTY 

5 ,  AA 
5.A8 
.7A 

0  00 
»*  •  W  b. 

2.67 
!  .73 
A.  AS 
32.  E3 


ANALYSIS: 

"52V26 

Vacuu.m  (P  k-‘a) 

=  AS.?8  (Deg  O 

•  SC. 00 
=  4  .  A  t  ( n '  s  ) 

=  3.8n8E+05  (W/m  2> 


PERCENT  UNCERTAINTY 

1  .AA 
'  .  57 
0.50 
0  .‘i3 
2.67 
.  66 
1  .A9 
2.31 


aiable  to  DTIC  does  nrH 

w  iegible  xeproducUon 
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D-UA  FOR  THE  UNCERTfilNTY  AMHLYSrS: 


Fjie  Harne: 

Pressure  Condition: 
Ctean  Temperature 
Uater  Fiou  Rate  ('/> 
Nater  Velocity 
Heat  Flux 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  r lou  Rate.  Md 
Reyri<-tids  Humber,  Re 
Heat  Flux.  Q 
._oci-Mean- T?n  DifF,  ._MTD 
kail  Resistance,  Ru 
Cve  rai  1  ri  . .  C  .  .  Uo 
Water-Side  H.T.C.,  Hi 
Stean-Sioe  rt.T.C..  Ho 


•:aTA  -nR  iHE  •■.•NCERTAINTY 

F-le  Vane: 

Pressure  Condition: 

Steam  Tomr-erature 
water  Rar^  <*/) 

water  VaicCiCy 
Heat  Flux 


UHCERTArrY  ANALYSIS: 
VARIABLE 

Mass  Fiou  Rate,  Md 
Reynolds  Humber,  Re 
Hear  Flux ,  q 
Locj-Mean- Ten  Diff.  LMTD 
Wall  Resistance.  Ru 
O'/erai  1  H  .  T  ,  C  .  .  Uo 
Water-Side  H.T.C..  Hi 
Steam-Side  H.T.C..  Ho 


F32AS8 

Atnospheric  < ! 0  I  kPa) 
*  3S  .  ‘33  ^  [>eo  C ) 

=  20,00 

» . IS  (n/s> 

-  6.502E  +  05  (W/ri’2) 


‘^FRCErr  UNCERTAINTY 

S  .AS 
S .  5C 
1  .7A 
0.07 


*>Cj  ^ 


ANALYSIS: 


Atnoapner ic  MQl  >  ^a> 
=  39.33  (Deg  O 

•  SO.  00 
=  A  .  :ir)  (  n  ) 

=  1.164r+06  (W/m'2) 


•■’ERCENT  UNCERTAINTY 

I  .AA 
I  .58 
G.A8 
0.  ]A 
2.S7 
SO 
1  .50 
2.80 
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Copy  available  to  DTIC  does  not 
peimit  fully  lepible  lepioductioa 


DATA  -V.U  THE  ijNCERTAINT /  hNALYSIS: 


E i 1  &  Name : 

F' reassure*  Condition: 
Sr-an  Tenperature 
Water  Flou  Rate  ('/.t 
Water  'Velocity 
Heat  Flux 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Nass  Flou  Rate.  Md 
Reynold?-  Number,  Re 
heat  Flax,  q 
Lcg-Ncan-”en  Diff.  LNTD 
Wail  Resistance.  Ru 
Overall  H.T.C..  Uo 
WatevSide  H.T.C..  Hi 
Steam“Side  H.T.C,.  Ho 


F39U33 

U  ac  uun  M  •  k  P  a  > 

•  48  50  < Dea  C) 

*  20.00 

=  !  03  <  m / ? ) 

-  2.579F.*05  (W/m'2) 


PERCENT  uncertainty 

5.42 
5.^5 
1.73 
f  .23 
4.39 
!  .75 

4.43 
13. C5 


DATA  “OR  The  uncertainty  ANALYSIS: 


File  '  /ane : 

“■  re? sure  Condition: 
3  team  Tenperature 
Hater  Flou  Rate  (71 
water  Vejocirv 
Heat  Flux 


FDSVSO 

Vacuum  (!’  fPa) 

=  43,4!  iDeq  C) 

-  A 0.00 

“  3.92  <  n/s  ' 

*  4,440E'*'0ti  (W/m  2> 


iiNCFRTAI.NTY  ANALYSIS: 


VARIABLE 

Mas?  Flou  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux,  q 
Log-Mean- Ten  Oiff.  LliTD 
Wall  Resistance.  Ru 
Overall  H  .  T  .  C .  ,  Uo 
Water-Side  H.T.C,.  Hi 
Steam-Side  H.T.C..  Ho 


PERCENT  UNCERTAINTY 

1.43 
I  .  5E) 

0.50 

0.48 

4.39 

.b? 

1.47 
1  .69 


<„o-aabU  10  DTIC  do.,  not 
11  InllY 


1  an 


FDP  THE  UUCERTHlMrY  HIJHLY5IS: 


File  fJame: 

^’ressare  Condition: 
Stoan  Temperature 
Water  Fiou  Rate  </l> 
Water  'Velocity 
heat  Flux 


UNCERTfilNTY  ANAlYSIS; 
v;ARIABLE 

Mass  F  lotu  Rate.  Md 
F^eynoids  Number.  Re 
heat  Flux,  q 
Log“Mean~ ^em  Dif^.  LMTD 
Wall  Resistance.  R'H 
Overall  H.T.C..  Uo 
Water “Side  H.T.C..  Hi 
Steam-Side  H.T.C. .  Ho 


DATA  FOR  the  UNCERTAINTY 

File  Name : 

Pressure  Condition: 

Steam  '•‘♦-mperature 
Water  ^ioi.i  Rate  <V. ' 
wa  r  ?  r  i oc . t  y 
Heat  F i ux 


LiNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  NiumDer,  Re 
Heat  Flux,  q 
Log-Mean-Iem  Diff.  LliTD 
Wall  Resistance,  Rw 
O'/eral  1  H.T.C.,  Uo 
Water-Side  H.T.C.,  Hi 
Steam-Side  H.T.C.,  Ho 


FCArtSe 

AtnospHeric  (!0*  '^Pa) 
=  A9.9f  (Deg  C) 

=  20.00 

*■  :  ,  OA  ( m/s  > 

“  8.131F+05  <W/m' 2) 


•^FRCEM"  IJNCEpTf^irnY 

5.2S 
5.  A 3 
1  .  72 
J.OR 
A . 

I  *  f  w 

A .  A2 
*8.0? 


ANALYSIS: 


~3‘3A9S 

Atmospneric  <101  kPa) 
99.73  (Deg  O 
SO  .00 
T.3A  (^(/'•.^) 

1.A09E+0?;  (W/n  2) 


'PERCENT  LINCCRTAINTY 


I  .A2 
1  .5A 
0.A8 
0.  !5 
A .  39 
.50 
1  .  Ab 
1  .89 
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